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MECHANISMS  OF  STRATOSPHERIC  OZONE  TRANSPORT 


Chapter  1 

INTRODUCTION  TO  THE  OZONE  TRANSPORT  PROBLEM 

1.1  Introduction 

Due  to  the  complex  nature  of  ozone  photochemistry  and 
atmospheric  dynamics  it  has  been  difficult  to  discern  the  mechanisms 
most  Important  for  the  maintenance  of  the  ozone  distribution.  To 
simultaneously  model  both  the  photochemistry  and  the  dynamics  is  a 
prodigious  task  as  each  system  involves  a  variety  of  temporal  and 
spatial  scales.  One  goal  of  atmospheric  scientists  is  to  understand 
transport  dynamics  so  that  they  may  be  simplified  and  then 
meaningful  predictions  of  ozone  perturbations  by  pollutants  can  be 
generated  for  a  considerably  reduced  computational  cost. 

The  oxygen  photochemistry  (Chapman,  1930;  Craig,  1965)  predicts 
that  the  maximum  amount  of  ozone  should  be  observed  in  the  tropical 
stratosphere  and  the  minimum  amount  observed  in  the  polar 
stratosphere.  However,  even  as  early  as  1930  it  was  known  that  the 
maximum  amount  of  ozone  occurs  at  extratropical  latitudes  (Dobson, 
1930),  and  it  thus  was  postulated  that  atmospheric  transport 
redistributes  ozone  from  its  tropical  source  regions  to  polar 
latitudes. 

Manuscript  approved  October  7, 1982. 


Simple  mean  cell  circulations  were  proposed  by  Brewer  (1949) 
and  Dobson  (1956)  which  were  capable  of  transporting  ozone  from 
tropical  to  polar  regions.  Planetary  wave  mechanisms  were  also 
suggested  that  were  capable  of  producing  poleward,  downward 
transport  (Newell,  1961,  1963).  As  more  data  became  available  and 
dynamical  theory  developed,  it  became  obvious  that  the  transport  of 
minor  constituents  was  a  complex  Interaction  of  planetary  waves  and 
the  mean  circulation  (Mahlman,  1969).  An  excellent  review  of  the 
history  of  stratospheric  circulation  and  transport  studies  is  given 
in  Hsu  (1980). 

Since  the  classic  study  of  Craig  (1950),  the  photochemical 
models  of  ozone  have  become  progressively  more  complex.  The 
remeasurement  of  rate  coefficients  by  Benson  and  Axworthy  (1957) 
indicated  that  oxygen-only  photochemical  models  would  predict  more 
ozone  than  is  observed  (Hunt,  1966).  Catalytic  destruction  cycles 
involving  hydrogen  oxides,  nitrogen  oxides,  and  halogen  oxides  have 
been  added  to  the  photochemistry  models.  With  the  addition  of  these 
reactions  it  became  obvious  that  anthropogenic  pollutants  might 
cause  serious  degradation  of  the  ozone  layer  (WMO,  1981).  Complete 
photochemical  models  now  require  the  use  of  over  100  chemical 
reactions  (Johnson  and  Fodolske,  1978).  Use  of  the  complete  set  of 
photochemical  equations  requires  that  the  dynamics  be  simplified  in 
order  to  make  the  problem  tractable.  The  most  common 
parameterization  has  been  to  treat  atmospheric  transport  as  if  it 


were  diffusive 


1.2  Diffusion  models 


The  dynamics  in  photochemistry  models  have  been  parameterized 
with  both  one -dimensional  (1-D)  and  two-dimensional  (2-D)  schemes. 
The  1-D  models  represent  transport  by  vertical  diffusion  from  high 
altitude  source  regions  to  the  lower  stratosphere  where  ozone  Is 
observed  to  accumulate  (McConnell  and  McElroy,  1973).  Hunten  (1975) 
pointed  out  that  vertical  diffusion  models  are  useful  because  of  the 
tendency  of  the  atmosphere  to  be  horizontally  stratified.  Upward 
motions  are  always  associated  with  compensating  downward  motions; 
therefore,  the  atmosphere  "looks  like  just  another  big  eddy" 
(Hunten,  1975).  Strobel  (1981),  however,  has  shown  that  1-D 
diffusion  models  are  only  strictly  valid  when  the  horizontal 
gradient  of  the  constituent  is  small,  a  condition  not  valid  for 
ozone . 

Two-dimensional  models  attempt  to  model  the  atmosphere  more 
realistically  by  approximating  both  vertical  and  latitudinal 
transport  (Hidalgo  and  Crutzen,  1977;  Harwood  and  Pyle,  1977;  Miller 
et  al . ,  1981).  In  2-D  models  a  mean  advective  field  is  usually 
prescribed  and  wave  motions  are  modeled  by  a  diffusion  scheme.  Such 
diffusion  parameterizatlons  are  often  based  on  the  empirical  theory 
of  Reed  and  German  (1965)  which  has  been  shown  to  be  inconsistent 
with  the  hydrodynamics  of  planetary  wave  transport  (Matsuno,  1980). 

Hidalgo  and  Crutzen  (1977)  and  Harwood  and  Pyle  (1977)  use  an 
Eulerian  mean  formulation  where  ozone  is  advected  with  a  zonal  mean 
velocity  and  the  zonal  mean  effects  of  the  planetary  waves  are 
modeled  by  diffusion.  In  another  2-D  model  Miller  et  al .  (1981)  use 
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Che  concepts  of  Che  Lagrangian-mean  circulation,  as  discussed  by 
Dunkerton  (1978),  Co  derive  a  Lagrangian-mean  circulation  for  the 
stratosphere.  Miller  et  al .  use  this  Lagrangian-mean  field  instead 
of  the  Eulerian  velocity  for  the  mean  field  advectlon. 

The  problems  of  applying  diffusion  models  to  atmospheric 
dynamics  has  been  discussed  by  Mahlman  (1975),  Tuck  (1979),  and  Pyle 
and  Rogers  (1980a).  These  authors  have  pointed  out  the  inadequacies 
of  diffusion  parameterizatlons  when  trying  to  model  countergradient 
transport,  the  spatial  variation  of  chemical  sources  and  sinks,  and 
the  diffusion  of  different  chemical  species.  These  processes  are 
important  for  ozone  in  some  parts  of  the  atmosphere  and  should  be 
properly  represented  in  ozone  models. 

Matsuno  (1980),  Clark  and  Rogers  (1978),  Dunkerton  (1980),  and 
Plumb  (1979)  have  shown  that  the  nature  of  planetary  waves  is 
advective  rather  than  diffusive  for  an  inert  tracer.  This  means 
that  the  diffusion  tensor  should  be  antisymmetric  rather  than 
symmetric  as  is  generally  assumed.  Matsuno  (1980),  Strobel  (1981), 
and  Holton  (1981)  show  that  there  are  also  diagonal  components  in 
the  diffusion  tensor  which  are  associated  with  chemical  source  terms 
and  wave  growth  and  decay. 

The  manner  in  which  diffusion  parameterizatlons  have  been 
applied  to  the  atmosphere  is  a  source  of  uncertainty  in  the  results 
of  chemical  models.  Pyle  and  Rogers  (1980a),  however,  argue  that  a 
dynamically  consistent  reformulation  of  the  diffusion  coefficients 
might  allow  for  the  continued  use  of  diffusion  models.  Strobel 
(1981)  has  formulated  such  a  consistent  theory  for  linear  waves  with 
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complex  frequencies.  This  theory  demonstrates  the  strong  coupling 


between  the  magnitude  of  the  diffusion  coefficients  and 

photochemistry.  Because  of  this  strong  coupling,  a  simple  diffusion 
model  can  adequately  predict  ozone  dynamics  only  under  special 
circumstances . 

Even  with  all  the  shortcomings,  diffusion  can  improve  the 
results  of  chemical  models.  A  recent  study  by  Pyle  and  Rogers 
(1980b)  shows  that  the  addition  of  diffusion  generates  more 

realistic  ozone  distributions  than  those  produced  by  the  diabatic 
circulation  alone.  Since  it  is  difficult  to  fully  model  both 
photochemistry  and  dynamics,  such  parameterizations  are  desirable 
and  necessary.  In  the  next  section  attempts  to  study  ozone  with  an 
accurate  representation  of  the  dynamics  and  a  simplified  chemistry 
will  be  discussed. 

1.3  General  circulation  models 

Several  general  circulation  models  have  been  used  to  study 
ozone  transport  and  the  interaction  of  radiation,  dynamics,  and 
photochemistry  (Clark,  1970;  Cunnold  et  al . ,  1975;  Schlesinger  and 
Mintz,  1979;  and  Mahlman  et  al.,  1980,  and  the  references 
therein).  These  models  emphasize  an  accurate  representation  of  the 
dynamics,  but  the  chemistry  is  modeled  in  a  simplified  manner. 
General  circulation  studies  have  met  with  some  success  in  modeling 
the  large  scale  features  of  the  ozone  distribution,  but  no  general 
circulation  model  has  ever  reproduced  all  of  the  large  scale 
features  of  the  ozone  distribution. 
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For  example,  the  quasigeostrophic  model  of  Clark  (1970) 
established  the  usefulness  of  the  Chapman  photochemistry  in  modeling 
ozone  and  demonstrated  that  the  stratospheric  circulation  was  able 
to  perturb  the  radiative  equilibrium  ozone  distribution  to  a 
distribution  that  approximated  observations.  In  particular, 

planetary  waves  were  found  to  transport  large  amounts  of  ozone  from 
equatorial  source  regions  to  high  latitudes.  By  increasing  the 
horizontal  gradient  of  temperature,  Clark  was  able  to  generate  high 
latitude  warmings;  however,  these  model  warmings  do  not  accurately 
represent  the  observed  characteristics  of  stratospheric  warmings. 
Therefore,  while  the  dynamics  in  this  model  do  indeed  cause  an 
increase  in  ozone  at  high  latitudes,  the  mechanisms  that  lead  to 
this  increase  are  in  doubt  due  to  the  inability  of  the  model  to 
accurately  represent  observed  dynamical  features. 

Using  the  same  dynamical  equations  of  Clark  (1970),  Cunnold  et 
al .  (1975)  generated  more  accurate  zonal  mean  circulations  than 

Clark;  however,  the  model  did  not  generate  stratospheric  warmings. 
In  their  formulation  of  the  ozone  continuity  equation  Cunnold  et  al . 
used  a  large  vertical  eddy  diffusion  coefficient,  based  on 
diffusivities  of  1-D  chemistry  models,  which  dominated  the  vertical 
transport  in  the  model.  Although  Cunnold  et  al .  did  calculate  a 
polar  spring  ozone  maximum,  the  accumulation  of  ozone  at  high 
latitudes  was  underestimated. 

Mahlman  et  al .  (1980)  used  a  3-D  primitive  equation  model  to 
perform  two  ozone  transport  experiments.  In  the  first  experiment 
the  ozone  was  specified  at  10  mb  (the  upper  boundary)  and  was  not 


allowed  to  deviate  from  this  climatological  value.  In  the  second 
experiment,  ozone  was  calculated  at  10  mb  using  a  simple  oxygen, 
hydrogen,  nitrogen  photochemistry.  Mahlman  et  al .  predicted  the 
vertical  structure  in  the  lower  stratosphere  quite  accurately  using 
either  of  the  specified  photochemistries.  This  indicates  that  the 
ozone  distribution  in  the  lower  conservative  region  of  the 
stratosphere  is  to  a  large  degree  independent  of  the  details  of  the 
photochemistry. 

The  most  serious  discrepancy  encountered  by  Mahlman  et  al .  was 
in  the  calculation  of  the  total  ozone  field.  Their  spring  maximum 
was  too  large,  occurred  too  far  south,  and  too  late  in  the  spring. 
The  model  of  Mahlman  et  al.  (1980)  did  not  produce  midwinter 
stratospheric  warmings  and,  indeed,  none  of  the  general  circulation 
models  accurately  represent  warming  events.  It  was  suggested  in 
Mahlman  et  al .  (1980)  that  the  transport  during  a  stratospheric 
warming  might  be  Important  in  improving  their  results. 

General  circulation  studies  are  so  complex  that  transport 
mechanisms  are  difficult  to  isolate  and  describe.  The  dynamics  are 
a  combination  of  diabatically  forced  circulations  and  wave  regimes, 
and  it  is  difficult  to  pinpoint  the  dominant  transport  mechanisms. 
Therefore,  the  use  of  simple  models  to  study  the  transport 
associated  with  specific  dynamic  mechanisms  might  be  useful  not  only 
for  application  to  chemical  models  but  for  interpreting  general 
circulation  model  results. 
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1.4  Lagrangian-mean  models 

In  general  Che  constituent  and  dynamics  fields  In  atmospheric 
studies  are  formulated  as  a  mean  field  and  a  deviation  field.  The 
mean  fields  are  generally  formed  by  taking  averages  around  a 
latitude  circle  in  what  is  called  the  Eulerian  formulation.  An 
alternative  to  the  Eulerian  zonal  mean  and  deviation  method  is  the 
Lagrangian-mean  method  (Andrews  and  McIntyre,  1978;  McIntyre,  1980), 
in  which  averages  are  taken  over  particular  fluid  tubes.  As  the 
averaging  operator  is  defined,  the  Lagrangian-mean  velocity  is  the 
center  of  mass  velocity  of  a  particular  ensemble  of  fluid  particles 
and  is  zero  when  the  conditions  of  the  noninteraction  theorem  (Boyd, 
1976;  Andrews  and  McIntyre,  1978,  Dunkerton,  1980)  are  met. 

With  the  above  attributes  it  may  seem  natural  to  study  tracer 
transport  in  the  Lagrangian-mean  system.  However,  it  turns  out  the 
Lagrangian-mean  transform  is  complicated  and  is  easily  formulated 
only  for  the  most  simple  flows  (McIntyre,  1980).  Furthermore, 
dispersion  of  the  particles  about  their  center  of  mass  may  render 
the  location  of  the  center  of  mass  of  little  use  in  the  actual 
location  of  the  tracer  (Hsu,  1980).  McIntyre  (1980)  has  suggested 
that  some  sort  of  continuous  reinitialization  of  the  lagrangian-mean 
model  might  solve  this  problem,  but  no  practical  scheme  is  currently 
available . 

Much  of  the  work  concerning  the  application  of  the  Lagrangian- 
mean  to  the  transport  of  tracers  in  the  stratosphere  is  based  on  the 
concepts  of  Dunkerton  (1978).  Dunkerton  shows  that  for  steady,  non¬ 
dissipating  planetary  waves,  the  Lagrangian-mean  circulation  of  the 
stratosphere  is  to  first  order  the  zonal  mean  circulation  first 
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estimated  by  Murgatroyd  and  Singleton  (1961).  This  circulation  is 
specifically  the  radiatlvely  forced  circulation  (diabatic 
circulation) . 

Planetary  waves  are  dissipative  and  not  steady  and  have  a 
significant  effect  on  tracer  transport;  therefore,  the  diabatic 
circulation  is  not  sufficient  to  describe  all  the  large  scale 
transport  processes.  Schoeberl  (1981a)  has  calculated  the 
Lagrangian-mean  circulation  associated  with  steady,  dissipative 
planetary  waves.  Schoeberl  found  that  planetary  waves  transport 
ozone  from  southern  to  northern  regions,  and  that  strong  downward 
transport  can  be  obtained  with  ‘‘realistic"  mean  zonal  wind  profiles. 
It  is  found  that  in  the  lower  stratosphere  the  Lagrangian-mean  flow 
associated  with  planetary  waves  is  dominate  over  the  diabatic 
circulation. 

The  2-b  model  of  Miller  et  al .  (1981)  used  a  modified  version 
of  Dunkerton's  diabatic  circulation  and  a  diffusion  tensor  based  on 
the  coefficients  of  Luther  (1974)  and  Hunten  (1975).  Since  the 
Lagrangian-mean  circulation  as  calculated  by  Dunkerton  (1978)  does 
not  include  the  effects  of  dissipative  planetary  waves,  the  use  of 
the  Luther  and  Hunten  diffusion  coefficients  does  not  accurately 
represent  planetary  wave  effects.  The  large  scale  transport  by 
waves  in  this  mode  1  are  advective,  and  any  diffusion 
parameterization  must  include  an  antisymmetric  component. 

In  another  approach  to  the  diffusion  problem  Holton  (1981)  and 
Strobel  (1981)  use  the  residual  circulation  defined  by  Andrews  and 
McIntyre  (1976)  to  calculate  tracer  transport.  Holton's  model 
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incorporates  both  planetary  wave  and  diabatic  effects.  A  diffusion 
model  is  then  used  to  calculate  the  small  scale  correction  to  the 
advection  by  the  residual  circulation.  This  formulation  involves  an 
antisymmetric  tensor  that  represents  an  advective  correction  to  the 
residual  circulation  which  is  assumed  to  be  negligible.  If  this 
antisymmetric  component  is  indeed  small  then  the  diffusive  part  of 
the  model  should  truly  represent  subplanetary  wave  effects. 

1.5  Objectives 

The  global  scale  transport  of  ozone  can  be  divided  into  two 
circulation  regimes:  transport  by  the  diabatic  circulation  and 
transport  by  large  scale  waves.  Each  of  these  transport  regimes  is 
important  and  both  must  be  modeled  in  a  complete  study.  There  are, 
of  course,  important  smaller  scale  processes  that  significantly 
affect  the  ozone  distribution.  These  small  scale  processes  are 
particularly  important  in  the  lower  stratosphere  and  in  the 
troposphere.  No  attempt  will  be  made  here  to  study  these  small 
scale  processes. 

In  this  dissertation  steady  state  and  time  dependent  planetary 
wave  models  are  used  to  investigate  ozone  transport  by  the  diabatic 
circulation  and  the  planetary  waves.  The  steady  state  model  is 
formulated  in  both  the  Eulerian  and  Lagrangian-mean  systems  and 
results  from  the  two  formulations  compared.  The  Lagrangian-mean 
circulation  is  also  compared  to  the  residual  circulation.  Using  the 
Lagrangian-mean  displacement  fields  a  Lagrangian  parcel  model  is 
developed  to  study  the  interaction  of  photochemistry  and  dynamics. 
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The  simplicity  of  the  steady  state  model  makes  it  easier  to  use  than 
the  time  dependent  model  to  understand  the  fundamental  transport 
mechanisms;  however,  the  steady  state  model  has  the  disadvantage 
that  potentially  important  transient  effects  are  neglected. 

In  the  time  dependent  model  quasigeostrophic  waves  are 
calculated  on  a  midlatitude  fi-plane .  This  model  is  the  same  as  used 
in  Schoeberl  (1982a,  1982b)  to  study  vacillation  and  wave 
saturation.  A  diabatic  circulation  based  on  the  flow  calculated  by 
Dunkerton  (1978)  is  superimposed  on  the  planetary  wave  circulation. 
With  this  model  it  is  possible  to  identify  the  transport 
characteristics  of  the  planetary  waves  and  the  diabatic  circulation 
and,  therefore,  measure  their  relative  importance. 

The  use  of  these  simple  models  aids  in  the  understanding  of  the 
specific  mechanisms  involved  in  the  interaction  of  photochemistry 
and  dynamics.  With  this  understanding  It  Is  possible  to  make  a 
rational  approach  to  the  development  of  dynamic  simplifications 
appropriate  for  chemical  models. 

In  Chapter  2  the  dynamic  and  transport  models  are  described 
with  particular  attention  given  to  the  transport  scheme.  In  Chapter 
3  the  transport  by  steady,  dissipative  waves  is  covered  and  the 
Eulerian,  Lagrangian-mean ,  and  residual  circulations  are  compared. 
The  transport  by  the  diabatic  circulation,  the  planetary  wave 
circulation,  and  the  two  circulations  combined  is  presented  for  the 
time  dependent  8-plane  model  in  Chapter  4.  Finally,  the  results  of 
the  various  models  are  discussed  and  interpreted  in  Chapter  5. 
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Chapter  2 


MODEL  DESCRIPTIONS  AND  THE  OBSERVED  OZONE  DISTRIBUTION 

2.1  Transport  model 

A  general  transport  model  has  been  developed  that  is  capable  of 
Interfacing  with  a  variety  of  dynamic  models  or  a  specified  time 
dependent  data  set.  The  transport  model  runs  in  an  off-line 
fashion;  that  Is,  the  constituent  distribution  does  not  affect  the 
dynamics.  Since  ozone  Is  radlatively  active,  with  the  use  of  an 
off-line  model  it  has  been  tacitly  assumed  that  the  change  in 
radiative  forcing  associated  with  the  changes  in  the  ozone  field  is 
small.  This  assumption  is  valid  as  most  of  the  transport  in  this 
model  takes  place  during  winter  when  insolation  is  small,  and  in  the 
infrared,  ozone  cooling  is  much  smaller  than  carbon  dioxide  cooling 
(Apruzese  et  al.,  1982).  Mahlman  and  Moxim  (1978)  discuss  in  detail 
the  differences  between  on-line  and  off-line  models. 

For  the  transport  model,  the  constituent  continuity  equation 
(Eq.  3.1)  is  written  in  discrete  form  using  fourth  order  centered 
finite  differences  in  the  meridional  and  vertical  directions.  The 
longitudinal  dependence  is  represented  by  zonal  harmonics,  and  the 
time  derivatives  are  estimated  with  second  order  centered 
differences  (leapfrog).  The  spatial  differences  are  written  in  the 
general  form 
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(2.1) 
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where  *  2/3  and  C£  “  -1/12  for  fourth  order  accuracy  (Ames,  1977; 
Haltiner  and  Williams,  1980).  Further  details  of  the  numerical 
scheme  are  given  in  Appendix  B. 

Mahlman  and  Sinclair  (1977)  review  in  some  detail  the 
difficulties  of  modeling  fluid  advection.  They  tested  several 
transport  algorithms  using  a  one-dimensional  model  with  a  steady, 
uniform  flow  and  cyclic  boundary  conditions.  The  most  serious 
deficiency  encountered  by  Mahlman  and  Sinclair  was  the  generation  of 
nonphysical  negative  mixing  ratios  which  had  to  be  removed.  This 
removal  process  is  called  "filling"  (Mahlman  and  Moxim,  1978). 
Filling  is  a  subjective  process  that  usually  involves  borrowing  mass 
from  surrounding  grid  points  in  a  mass  conserving  fashion  so  that 
negative  mixing  ratios  are  increased  to  zero. 

Mahlman  and  Sinclair  found  that  the  pseudo-spectral  method 
produced  the  best  results  when  a  large  number  of  waves  and  a  short 
time  step  were  used.  However,  the  pseudo-spectral  scheme  requires  a 
large  amount  of  computer  time  and  memory.  Of  the  other  schemes 
tested  the  fourth  order  derivative  scheme  worked  well  and  is 
computationally  efficient  (see  also  Kreiss  and  Oliger,  1972). 

Clancy  (1981)  has  suggested  that  the  advection-dif fusion 
equation  can  be  efficiently  integrated  for  a  wide  range  of 
parameters  using  second  order  spatial  differences  and  a  forward  time 
scheme.  A  comparison  of  the  current  model  with  Clancy's  scheme 
illustrates  the  computational  gain  of  the  fourth  order  differences. 

13 


Following  Clancy,  the  advection-diffusion  equation, 


3S 

at 


as  „  a2s 

"U  ^  +  D  —  2 
3x 


is  integrated  with  At  =  3600  sec.  Ax  =  10^  ra,  u  *  .6  m  sec-*,  and  D 
«  10^  m^  sec-*.  The  initial  condition  is 


S  =  SQsin(kx)  (t  -  0) 

and  the  boundary  conditions  are 

S  =  -SQexp(-k^Dt )sin(kx-ukt ) 

for  x  ■  0  and  10*  m,  where  k  =  2ir ( 10  m)  .  The  results  after  30 
days  using  second  order  spatial  differences  and  forward  time 
differences  are  shown  in  Fig.  2.1a  (same  as  Clancy  Fig.  2).  In  Fig. 
2.1b  forward  time  differences  are  still  used,  but  the  second  order 
differences  are  replaced  with  the  fourth  order  estimates.  The  phase 
error  in  Fig.  2.1b  has  been  reduced  to  nearly  zero,  but  as  in  Fig. 
2.1a  the  amplitude  is  overestimated. 
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(same  as  Clancy  (1981),  Fig.  2) 


c.  Centered  In  time,  fourth  order  in  space 

Figure  2.1.  A  comparison  of  the  analytic  and  numerical  solution 
of  the  advection-dif fusion  equation  using  different  numerical  schemes. 

(  +  on  the  numerical  solution) 
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In  Fig  2.1c  the  same  problem  has  been  redone  using  centered 
time  differences  and  fourth  order  spatial  differences.  Both  the 
phase  error  and  the  amplitude  error  have  been  reduced  to  the  point 
that  the  numerical  solution  and  the  analytic  solution  are  virtually 
indistinguishable.  The  use  of  forward  time  differences,  as  in  Fig. 
2-la  and  2.1b,  leads  to  an  unstable  solution  when  the  diffusion  is 
zero  (Haltiner  and  Williams,  1980,  p.  130).  The  incorporation  of 
diffusion  by  Clancy  serves  to  stabilize  his  scheme.  A  more  detailed 
analysis  of  this  problem  is  presented  in  Appendix  B. 

Using  second  order  differences  without  diffusion  to  model  ozone 
transport  leads  to  unrealistic  results  because  2Ax  waves  generated 
in  regions  of  large  gradients  propagate  throughout  the  domain  and 
dominate  the  computed  transport .  Only  with  the  incorporation  of 
very  large  diffusion  can  this  effect  be  countered  for  long 
integrations.  The  inferiority  of  second  order  transport  models  is 
well  illustrated  in  Mahlman  and  Sinclair  (1977). 

Besides  filling,  there  are  other  nonphysical  processes  that 
must  often  be  incorporated  into  transport  models  in  order  to  obtain 
reasonable  numerical  results.  Diffusion  is  commonly  introduced  to 
represent  subscale  effects  (Mahlman  and  Moxlm,  1978).  As  indicated 
above,  diffusion  generally  increases  the  numerical  stability  of  a 
finite  difference  scheme  and  reduces  the  occurrence  of  negative 
mixing  ratios. 

At  the  boundary  some  technique  must  be  used  in  order  to  close 
the  finite  difference  scheme.  In  this  model  either  the  ozone 
density  or  flux  is  specified  at  the  boundary,  and  the  derivative  at 


the  first  interior  point  is  represented  by  second  order  centered 
differences  (C^»  1/2,  C2“  0  in  Eq.  2.1).  With  this  closure  system 
the  boundary  points  and  the  first  two  interior  points  are  directly 
affected  by  the  prescription  of  the  boundary. 

The  closure  method  chosen  generally  affects  the  ability  of  the 
model  to  conserve  mass  due  to  the  uncompensated  difference  terms  at 
each  boundary.  In  a  problem  where  there  is  no  flux  through  the 
boundary  and  no  source  or  sink  terms,  it  is  possible  to  calculate 
the  mass  change  caused  by  these  uncompensating  terms  and  then 
distribute  an  appropriate  mass  correction  over  the  model  domain.  In 
a  problem  with  boundary  fluxes  and  sources,  mass  correction  becomes 
more  complex. 

Closure  at  the  boundaries  has  another  effect  which  particularly 
manifests  itself  in  the  ozone  problem.  Since  the  mixing  ratio 
estimates  at  or  near  the  boundaries  are  less  accurate  estimates  than 
at  the  interior  points,  transport  into  the  polar  region,  a  region  of 
special  interest,  will  be  Inaccurate. 

The  philosophy  taken  in  this  study  is  to  incorporate  as  little 
subjectivity  as  possible,  so  that  Important  transport  mechanisms 
will  be  revealed.  Neither  diffusion  nor  mass  correction  have  been 
incorporated  Into  this  model,  and  it  has  not  been  necessary  to 
Incorporate  filling. 

The  omission  of  diffusion  means  that  the  lower  regions  of  the 
stratosphere  and  troposphere  are  expected  to  be  poorly  modeled.  In 
particular  there  is  no  mechanism  for  tropospheric-stratospheric 
interaction  and  there  are  no  subplanetary  scale  dynamics  to  produce 
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small  scale  mixing  effects.  As  a  result,  there  Is  sometimes  a 
tendency  for  ozone  to  build  up  to  unrealistic  levels  in  the  lower 
stratosphere.  In  certain  instances  an  ozone  "shock  front”  develops 
and  the  Unite  difference  estimates  of  the  derivatives  become  very 
poor,  ultimately  leading  to  the  failure  of  the  model  in  all 
regions.  Strong  ozone  fronts  do  not  develop  with  the  dynamics  used 
here . 

The  ability  of  the  model  to  conserve  mass  has  been  studied 
using  an  inert  tracer.  The  model  does  not  conserve  mass  exactly  and 
the  entirety  of  the  mass  correction  can  be  attributed  to 
noncompensating  finite  difference  terms.  The  mass  discrepancy 
remains  small  for  the  integration  times  considered  here.  When  ozone 
chemistry  is  included  the  total  ozone  in  the  model  domain  never 
differs  unrealistically  from  the  initial  total  ozone.  For  these 
reasons  it  was  judged  that  the  incorporation  of  any  mass  correcting 
techniques  would  contaminate  the  model  results  at  least  as  badly  as 
the  uncorrected  mass.  Therefore,  the  results  presented  here  have  no 
mass  correction  added. 

2.2  Chemistry  model 

The  chemistry  used  in  this  study  is  extremely  simple,  namely: 

Q  -  -  Am  (A  >  0)  (2.2) 

where  Q  is  the  photochemical  source  term  (see  Eq.  3.1),  A  is  the 
photochemical  relaxation  rate,  and  m  is  the  constituent  mixing 


ratio.  This  formulation  requires  that  the  perturbations  relax  back 


to  zero,  while  changes  in  the  mean  field  relax  to  the  mean 
equilibrium.  The  rate  of  relaxation  is  inversely  proportional  to  X, 
which  is  assumed  to  be  a  function  of  height  only.  The  X  values  are 
based  on  those  of  Hartmann  and  Garcia  (1979). 

The  chemical  reactions  that  the  X  values  are  based  on  are  the 
oxygen-only  reactions  (Chapman  reactions)  with  a  modified  ozone 
destruction  rate  to  represent  catalytic  destruction.  Ozone 
chemistry  is  temperature  dependent;  however,  the  temperature 
dependence  is  weak  (Hartmann  and  Garcia,  1979;  Garcia  and  Hartmann, 
1980)  and  will  be  neglected  here. 

It  is  convenient  to  divide  the  atmosphere  into  three  regions  in 
order  to  discuss  ozone  transport  (Fig.  2.2).  The  lower  region  is 
the  region  in  which  ozone  is  assumed  to  be  a  conservative 
tracer  (X  «  ku).  In  the  upper  region  ozone  photochemistry  is  so 
fast  (X  »  ku  )  that  the  transport  is  unimportant.  The  middle 
region  is  defined  by  the  ratio  .1  <  X/uk  <  10  and  is  called  the 
transition  region  (Hartmann  and  Garcia,  1979).  The  Interaction  of 
the  photochemistry  with  the  dynamics  is  expected  to  be  greatest  in 


the  transition  region. 


Figure  2.2.  Schematic  representation  of  the  atmosphere  divided 
into  photochemical,  transitional,  and  conservative 
regions  f ased  on  the  ratio  of  the  dynamical  to  the 
photochemical  time  scale. 
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2.3  Dynamic  models 

The  planetary  waves  are  represented  in  all  the  dynamic  models 
by  the  quasigeostrophic  potential  vorticity  equation.  The  governing 
equation  is 


+  a  +  u 
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where  q'  is  defined  as  in  Holton  (1975).  The  mean  flow  equation  is 
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(2.4) 


q  is  the  wave  induced  change  in  the  background  potential  vorticity 
(Jq.  Damping  is  modeled  as  Newtonian  cooling  and  Rayleigh  friction, 
which  are  assumed  to  be  equal,  and  is  represented  by  a.  F  is 
diabatic  forcing  other  than  Newtonian  cooling.  Other  nomenclature 
is  given  in  the  Appendix  A. 


a.  Steady  state  model 

In  the  steady  state  model  is  assumed  to  be  zero  and  0-plane 

geometry  is  used.  The  damping  is  a  slowly  varying  function  of 
height  except  in  the  uppermost  regions  where  it  is  assumed  to 
increase  at  an  exponential  rate.  This  region  of  increase  provides  a 
sponge  layer  to  reduce  reflections  from  the  upper  boundary  and  is 
well  above  the  region  where  ozone  transport  is  being  studied.  The 
time  scale  of  the  damping  below  the  sponge  layer  is  approximately  10 
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days,  which  is  representative  of  the  physical  processes  in  the 
middle  or  upper  stratosphere  (Schoeberl  and  Strobel,  1978; 
Dickinson,  1973).  The  diabatic  forcing,  F,  is  assumed  to  be  zero; 
therefore,  only  the  steady  state  planetary  wave  transport  is 
modeled . 

Equations  (2.3)  and  (2.4)  are  solved  for  a  wavenumber  one 

perturbation  with  forcing  analogous  to  a  mean  zonal  wind  of  10  m 
sec~*  flowing  over  a  1  km  mountain.  The  mean  zonal  wind  profile  and 
the  geopotential  amplitude  are  shown  in  Fig.  2.3.  The  mean  zonal 
wind  is  presumed  to  approximate  the  polar  night  jet.  The 

geopotential  reaches  a  maximum  amplitude  of  about  1  km  in  the 

mesosphere  and  does  not  exceed  the  steady  state  saturation  limit 
(Schoeberl,  1982a).  The  magnitude  of  the  geopotential  is  in  rough 
agreement  with  observation  (van  Loon  et  al.,  19  73).  The  eddy 
geopotential  field  is  used  to  calculate  the  wave  forced  mean 

meridional  and  vertical  velocity. 

It  is  assumed  in  the  steady  state  model  that  wave  effects  on 
the  mean  flow  are  not  so  large  that  the  wave  field  needs  to  be 

recomputed  with  the  perturbed  mean  flow.  This  restricts  the 

validity  of  the  steady  state  model  to  small  amplitude  waves  and  is 
consistent  with  the  transport  approximations  to  be  made  later 

(Section  3.2).  This  model  is  the  Eulerian  analogue  to  the 

Lagranglan-mean  model  of  Schoeberl  (1981a). 
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b.  Magnitude  of  the  steady  state  geopotential. 
Figure  2.3. 


HEIGHT  (km) 
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U  (m  sec-1) 


for  the  time  dependent  f}-plane  model. 
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b.  Time  dependent  6-plane  model 

In  the  time  dependent  6-plane  model  the  diabatlc  forcing,  F  in 
Eqs.  (2.3)  and  (2.4),  Is  assumed  to  be  zero,  and  the  Newtonian 
cooling  and  Rayleigh  friction  are  similar  to  that  of  the  steady 
state  model.  Diabatlc  effects  are  modeled  by  adapting  a  circulation 
from  Dunkerton  (1978)  and  Murgatroyd  and  Singleton  (1961)  and 
superimposing  this  circulation  on  the  calculated  planetary  wave 
circulation. 

The  waves  are  forced  by  zonal  flow  over  a  mountain.  The  basic 
state  mean  zonal  wind  profile  is  shown  in  Fig.  2.4.  This  model  is 
quasilinear  in  that  after  every  time  step  the  wave  perturbed  mean 
field  is  calculated;  the  waves  are  then  recalculated  using  the 
updated  mean  fields.  The  mountain  height  is  a  model  parameter  and 
the  behavior  of  the  model  varies  depending  on  the  mountain  height. 
The  planetary  wave  model  is  described  in  Schoeberl  (1982b). 

For  small  mountains  equilibrium  is  reached  in  which  westerlies 
are  present  everywhere.  For  mountains  in  an  intermediate  range  the 
model  produces  at  least  one  stratospheric  warming  and  then  is 
observed  to  reach  a  steady  state  with  a  permanent  layer  of 
easterlies.  For  large  topography  the  model  produces  a  series  of 
warmings,  i.e.  it  vacillates.  For  the  transport  experiments  two 
mountain  heights  will  be  used:  600  m  and  900  ra.  Both  of  these 
mountain  heights  are  in  the  intermediate  range  where  there  is  a 
warming  and  the  eventual  onset  of  a  permanent  layer  of  easterlies. 

The  8-plane  transport  model  is  centered  at  60°  N  and  has  a 
horizontal  extent  of  5000  km.  This  corresponds  to  a  latitudinal 
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extent  from  37.5°  to  82.5°.  The  domain  extends  from  the  ground  to 
100  km.  There  are  20  grid  points  in  the  horizontal  and  40  grid 
points  in  the  vertical  and  a  perimeter  of  points  around  the  grid  at 
which  the  boundary  values  are  prescribed. 

At  the  top  and  bottom  of  the  model  the  constituent  flux  is  set 
equal  to  zero.  This  boundary  condition  is  accurate  at  the  top 
since  the  oz^ne  is  in  photochemical  equilibrium.  At  the  northern  and 
southern  boundaries  the  flux  associated  with  the  planetary  waves  is 
assumed  to  be  zero.  At  the  southern  boundary,  where  ozone  is 
observed  to  have  relatively  little  annual  variation,  the  flux  by  the 
diabatic  circulation  is  calculated  assuming  that  ozone  is 
constant.  The  flux  associated  with  the  diabatic  circulation  at  the 
northern  boundary  is  assumed  to  be  zero . 

The  closure  of  the  finite  difference  scheme  directly  involves 
the  prescribed  boundary  point  and  the  first  interior  point.  The 
second  and  third  interior  point  are  directly  associated  with  the 
first  interior  point  at  each  time  step.  Because  of  this  possible 
direct  contamination  by  the  closure  system,  results  will  be 
presented  at  horizontal  grid  points  4  (46°)  and  17  (74°)  as 

representative  of  the  northern  and  southern  regions  of  the  model . 
Table  2.1  lists  the  grid  points,  the  corresponding  distance  from  the 
southern  bondary  of  the  8-channel  ,  and  the  latitude  for  the  time 
dependent  model. 


Table  2.1 

Horizontal  Distance  Scales  for  the  Tine 
Dependent  0-plane  Model 


Grid  Point 

Y  (km) 

Latitude  0 

1 

235 

39.6 

2 

476 

41.8 

3 

714 

43.9 

4 

952 

46.1 

5 

1190 

48.2 

6 

1428 

50.4 

7 

1667 

52.5 

8 

1905 

54.6 

9 

2143 

56.8 

10 

2381 

58.9 

11 

2619 

61.1 

12 

28  5  7 

63.2 

13 

3095 

65.4 

14 

3333 

67.5 

15 

35  71 

69.6 

16 

3810 

71.8 

17 

4048 

73.9 

18 

4285 

76. 1 

19 

4523 

78.2 

20 

4  761 

80.4 

27 


2.4  Observations  of  ozone 

In  order  to  evaluate  the  success  of  the  transport  model  it  is 
necessary  to  compare  the  calculated  results  to  the  observed  ozone 
distribution.  In  this  section  a  brief  review  of  the  observed 
behavior  of  ozone  will  be  given.  Some  mention  should  be  made 
concerning  the  units  used  in  reporting  the  observations.  Hie 
transport  calculations  will  be  done  using  the  mixing  ratio  (ppm)  as 
the  basic  unit  of  measure .  Many  of  the  observations  have  been  made 
in  other  units  such  as  ozone  partial  pressure  or  absolute  density. 
Due  to  the  variation  of  atmospheric  density,  the  ozone  distribution 
appears  quite  different  in  the  various  units.  Kreuger  and  Minzner 
(1976)  present  conversion  factors  in  their  Table  2.  Much  of  the 
general  review  of  the  ozone  distribution  is  taken  from  Craig  (1965) 
and  Du'tsch  (19  74). 

The  ozone  distribution  in  the  Southern  Hemisphere  is  less  well 
known  and  considerably  different  from  that  in  the  Northern 
Hemisphere.  The  distribution  in  the  Southern  Hemisphere  will  not  be 
discussed  here  since  the  dynamics  are  more  appropriate  for  the 
Northern  Hemisphere.  There  will  be  some  comment  made  about  transport 
mechanisms  in  the  Southern  Hemisphere  based  on  model  results  in 
Chapter  5. 

At  extratropical  latitudes  the  ozone  density  reaches  a  maximum 
in  the  lower  stratosphere  between  15  and  30  km.  The  mixing  ratio 
maximum  occurs  at  somewhat  higher  altitudes  between  35  and  45  km. 
Above  and  below  the  ozone  maximum  the  average  state  shows  a 
monotonic  decrease  of  ozone  concentration;  though,  at  any  time  an 
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ozone  sounding  night  show  considerable  vertical  structure  due  to 
small  scale  transport.  There  is  ozone  of  stratospheric  origin  in 
the  troposphere,  but  the  amounts  are  quite  small;  and  there  is  no 
attempt  made  here  to  model  tropospheric  ozone. 

The  absolute  maximum  amounts  of  ozone  in  the  Northern 
Hemisphere  occur  at  polar  latitudes  in  the  lower  stratosphere  during 
late  winter  and  early  spring.  This  is  in  complete  disagreement  with 
photochemical  theory  which  predicts  maximum  amounts  in  the  tropics 
at  high  altitudes  in  the  middle  and  upper  stratosphere.  Figure  2.5 
shows  a  March/April  pole  to  pole  cross  section  of  ozone  in  both 
ozone  density  and  ozone/air  number  mixing  ratio.  The  predominant 
maximum  seen  at  high  latitudes  in  the  density  field  is  not  clearly 
seen  in  the  mixing  ratio  field  since  this  maximum  occurs  in  a  region 
of  relatively  high  atmospheric  density.  Two  important  features  to 
note  in  the  mixing  ratio  field  are  the  positive  horizontal  gradient 
at  high  latitudes  in  the  lowest  part  of  the  stratosphere  and  the 
negative  horizontal  gradient  at  high  latitudes  in  the  middle  and 
upper  stratosphere  (Dutsch,  1974). 

In  the  photochemical  region  the  variation  of  ozone  is  quite 
small.  In  the  lower  statosphere,  however,  the  variation  is  large  on 
both  daily  and  seasonal  time  scales.  Much  of  the  daily  variation  is 
due  to  the  infusion  of  synoptic  size  weather  systems  into  the 
stratosphere  and  will  not  be  modeled  here.  The  seasonal  variation 
is  due  to  the  diabatic  circulation  and  planetary  wave  transport  and 
will  be  investigated  in  this  dissertation. 
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Figure  2.6*  6  year  average  time-height  cross  section 

of  ozone  partial  pressure  at  Arosa,  Switzerland 
(from  I>Jtsch,  19  74). 
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a.  Arctic  type 

b.  Polar  type 

c-  Subtropical  type 
Temperature  on  left  in  °C. 

Figure  2.7  Vertical  distribution  of  ozone  from  3  soundings 
in  late  winter  and  early  spring  (ymb)  (from  Dutsch,  1974). 
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Figure  2.6  shows  a  time  cross  section  of  vertical  ozone  over 
Arosa,  Switzerland.  The  most  important  feature  to  note  is  the 
strong  increase  between  30  and  50  mb  (21-25  km)  which  begins  in 
December  and  reaches  a  maximum  in  late  winter.  By  the  first  of 
April  there  is  a  definite  decreasing  trend.  Since  this  maximum 
appears  in  the  conservative  portion  of  the  stratosphere,  it  can  be 
deduced  that  transport  is  responsible  for  the  ozone  buildup. 
Another  notable  feature  is  the  minimum  concentration  that  occurs  at 
the  level  of  the  tropopause . 

In  Figure  2.7  is  shown  the  vertical  distribution  from  three 
soundings  in  late  winter  and  early  spring.  In  the  polar  type 
sounding  (b)  there  is  a  pronounced  secondary  density  maximum  in  the 
lowest  part  of  the  stratosphere  and  a  primary  maximum  around  21 
km.  In  the  Arctic  type  (a)  there  is  just  one  broad  maximum  centered 
at  100  mb.  In  the  subtropical  type  profile  (c)  there  is  only  one 
large  maximum,  and  there  is  relatively  little  ozone  in  the  lower 
stratosphere. 

Retrieval  of  vertical  profiles  in  general  requires  the 
inversion  of  remote  sensing  observations  which  involve  strongly 
overlapping  contribution  functions.  Consequently,  there  are 
significant  inaccuracies  in  the  determination  of  the  vertical 
profile.  Ozonesondes  have  been  used  to  directly  determine  the 
vertical  structure,  but  these  observations  are  limited  in  both  space 
and  time.  Total  ozone  has  been  studied  longer  and  is  known  with 
greater  accuracy  than  the  vertical  structure;  therefore,  total  ozone 
will  be  used  in  many  of  the  comparisons  of  the  model  results  to 
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observation.  Using  total  ozone  for  comparison  has  the  disadvantage 
that  model  errors  might  be  masked  during  the  Integration  to  form  the 
total  amounts  - 

Many  of  the  features  of  the  behavior  of  total  ozone  have  been 
known  since  the  early  days  of  ozone  observations.  Using  a  network 
of  only  six  instruments,  Dobson  (1930)  deduced  that  both  total  ozone 
and  the  temporal  variation  of  total  ozone  increased  with  increasing 
latitude.  Perhaps  the  most  extensive  documentation  of  analyzed 
total  ozone  observations  is  the  atlas  of  London  et  al.  (1976)  which 
contains  ten  years  of  data  from  1957-1967.  The  average  latitude¬ 
time  cross  section  of  total  ozone  in  Dobson  Units  (DU)  is  shown  in 
Fig.  2-8  and  is  not  greatly  different  in  the  Northern  Hemisphere 
from  the  cross  section  published  by  Gotz  (1951).  This  figure  shows 
a  rapid  buildup  of  ozone  at  middle  and  high  latitudes  during 
winter.  The  maximum  is  reached  at  polar  latitudes  around  the  first 
of  March.  During  early  winter  the  maximum  amount  of  ozone  is  found 
to  between  60°  and  70°  N.  There  is  a  lack  of  data  north  of  80°  so 
that  the  results  in  the  extreme  north  are  of  dubious  quality. 

The  regularity  of  the  total  ozone  field  is  demonstrated  in  Fig. 
2.9  taken  from  Hilsenrath  et  al  ♦  (1979).  These  data  from  the  NIMBUS 
4  BUV  experiment  show  some  fluctuation  in  the  annual  variation  of 
total  ozone,  but  the  large  scale  features  remain  quite  constant 
during  the  time  series. 

Figure  2.10  from  Zullig  (1973)  further  demonstrates  the  extent 
of  the  Interannual  variation  in  the  total  oz one  field.  In  this 
figure  the  average  amount  of  total  ozone  in  the  Northern  Hemisphere 
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north  of  40°  for  three  winters  is  shown.  Even  with  the  large  amount 
of  averaging  involved  in  the  formation  of  these  graphs,  there  is  an 
obvious  temporal  wave  evident  in  the  data.  In  the  two  winters 
during  which  warmings  occur  (1962/63,  1967/68),  the  maximum  amount 
of  ozone  is  approximately  62  larger,  and  the  increase  is  more 
sudden,  than  in  the  winter  during  which  no  major  warming  occurs. 
Even  though  Zullig's  observations  are  a  rather  limited  amount  of 
data  on  which  to  base  transport  theories,  it  is  safe  to  conclude 
that  major  stratospheric  warmings  have  a  strong  influence  on  the 
ozone  distribution. 
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MONTH 

10  year  average  latitude-time  cross  section 
of  total  ozone  in  Dobson  Units  (DU) 

(from  London  et  al . ,  1976). 
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Figure  2.11  Reproduction  of  Figure 


from  mid-December  thro 
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The  basic  features  of  the  ozone  variabilty  that  will  be  studied 
here  are  the  total  ozone  amount  at  extratropical  latitudes  and  the 
transport  of  ozone  into  the  lower  stratosphere  during  winter.  Total 
ozone  increases  at  all  latitudes  during  winter  with  the  magnitude  of 
the  increase  being  largest  near  the  pole.  During  early  winter  the 
maximum  is  in  subpolar  latitudes.  Much  of  this  increase  in  total 
ozone  is  associated  with  the  injection  of  large  amounts  of  ozone 
between  21  and  25  km.  There  is  little  change  in  the  ozone  above  40 
km  and  the  observed  changes  in  this  region  can  be  explained 
photochemically .  The  total  ozone  reaches  a  maximum  in  early  March 
at  polar  latitudes  and  decreases  uniformly  until  it  reaches  a 
minimum  in  the  fall. 

The  time  dependent  0-plane  transport  model  will  be  used  to  run 
simulations  from  mid-December  until  April.  In  order  to  facilitate 
comparisons  to  observations  a  reproduction  of  the  appropriate 
portion  of  the  average  total  ozone  field  given  by  London  et  al . 
(1976)  is  presented  in  Fig.  2.11. 
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Chapter  3 


TRANSPORT  BY  STEADY  STATE  PLANETARY  WAVES 

3*1  Introduction 

In  this  chapter  steady  state  transport  by  stationary  planetary 
waves  on  a  fl-plane  is  presented-  In  Section  3.2  the  continuity 
equation  is  developed  for  small  amplitude  waves  in  both  the  Eulerian 
and  Lagrangian-mean  formulations,  and  the  relation  between  the 
residual  circulation  and  the  Lagrangian-mean  velocity  is  derived. 
In  Section  3.3  model  results  are  presented  and  a  Lagrangian  model 
based  on  particle  trajectories  is  given.  The  results  are  discussed 
and  summarized  in  Section  3.4. 

3.2  Constituent  dynamics 

The  ozone  continuity  equation  is  formulated  in  both  the 
Eulerian  and  Lagrangian-mean  systems.  In  general,  the  constituent 
continuity  equation  Is 


dm 

Dt 


-  Q 


(3.1) 


where  Q  represents  the  source  term.  By  decomposing  the  fields  into 
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zonal  averages  and  deviations  and  linearizing,  the  Eulerian 
formulation  is  obtained: 


in. 

3t 


+  u. 


=  Q 


(3.2) 


and 
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(3.3) 


where  repeated  indices  indicate  summation.  It  is  pertinent  to  note 
that  the  zonal  mean  equation  contains  eddy  flux  terms  which 
represent  the  effects  of  wave  transport.  The  tendency  for  the  eddy 
transport  to  be  counteracted  by  the  advection  by  the  mean  flow  has 
been  discussed  extensively,  and  represents  a  major  difficulty  of 
atmospheric  transport  problems  (Dunkerton,  1980;  Mahlman  and  Moxim, 
1978). 

Using  the  Lagrangian-mean  formulation  (Andrews  and  McIntyre, 
1978)  it  is  possible  to  redefine  the  averaging  operator  such  that 
wave  forcing  does  not  appear  in  the  mean  equation.  Using  this 
concept  the  mean  continuity  equation  is  written  as 


+  u 


L  3 
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(3.4) 


In  this  system  the  time  rate  of  change  of  the  Lagrangian-mean  ozone 
distribution  is  equal  to  the  advection  by  the  Lagrangian-mean 
velocity  plus  a  source  term. 
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The  difference  between  the  Lagrangian-mean  and  the  Eulerian 
zonal  mean  is  defined  as  the  Stokes  correction 

s  _  L 
M  =  U  ~  \i 

Relating  the  Lagrangian-mean  transform  to  the  Eulerian  zonal  mean  by 
the  premise  that  the  unperturbed  fluid  would  be  described  by  the 
Eulerian  flow,  it  is  possible  to  define  the  Stokes  correction  to 
second  order  as 
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(3.5) 


where  ^  *  (£lf  52»  S3)  ”  (S*  n,  C)  are  the  displacements  from  the 
equilibrium  position  of  the  fluid. 

The  displacement  fields  are  related  to  the  Eulerian 
perturbation  velocity  fields  by  the  following  equation  (Matsuno  and 
Nakamura,  1979): 


+  5i  -u;  +  <j|j  :i  +  —  (3-6> 

The  mean  meridional  circulation  is  forced  by  the  waves  and  is 
assumed  to  be  second  order  in  wave  amplitude  (0(a2)).  Therefore, 
the  changes  in  the  mean  ozone  field  are  O(a^)  and  the  terms 
involving  v,  w  in  (3.3)  and  (3.6)  can  be  ignored.  The  perturbation 
equation  can  now  be  written  to  O(a^)  as 
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Using  the  mass  continuity  equation,  (3.2)  may  be  written  as 
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where  y  now  represents  changes  in  the  mean  field  from  the 
equilibrium  value  yQ. 

The  relation  between  the  perturbation  velocity  fields  and  the 
displacement  fields  (Eq.  3.6)  can  be  rewritten  for  steady  waves  as 


ikilc  -  «•'  +  ng  + 


ikuri  ■  v" 
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These  are  the  equations  used  to  calculate  the  Lagrangian-mean 
displacements . 

Under  the  assumptions  above,  the  following  relationships 
between  the  Eulerlan  and  the  Lagrangian-mean  quantities  can  be 
written  for  steady  waves: 

+  O(a-)  (3.9s) 
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Therefore  (3.4)  can  be  written  to  the  same  accuracy  as  (3.8)  as 
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(3.10) 


The  first  term  on  the  right  hand  side  of  (3.10)  will  be  referred  to 
as  the  advective  term;  the  second  term  is  the  "stirring"  term;  and 
the  third  term  is  the  mean  chemistry.  The  stirring  term  represents 
the  effects  of  the  planetary  waves  moving  ozone  parcels  through  a 
region  of  varying  photochemistry .  This  term  will  be  discussed  in 
detail  later.  Equations  (3.7),  (3.8),  and  (3.10)  are  the  ozone 
continuity  equations  used  in  this  study. 

A  number  of  authors  have  suggested  that  the  Lagrangian-mean 
velocity  may  be  a  suitable  quantity  to  approximate  the  large  scale 
transport  of  ozone  in  chemical  models  (Dunkerton,  1978;  Miller  et 
al .,  1981;  McIntyre,  1980).  However,  since  the  Lagrangian-mean 

velocity  can  be  difficult  to  calculate,  it  has  also  been  suggested 
that  the  residual  circulation  defined  by  Andrews  and  McIntyre  (1976) 
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night  be  an  adequate  approximation  for  the  Lagranglan-mean  velocity 
(Holton,  1981).  Strobel  (1981)  has  derived  a  general  linearized  2-D 
transport  model  using  the  residual  circulation  for  waves  with 
complex  frequencies.  The  relation  between  the  residual  velocity  and 
the  lagrangian-mean  velocity  for  steady,  dissipating  planetary  waves 
is  outlined  below. 

The  eddy  mass  continuity  equation  can  be  written  in  terms  of 
the  parcel  displacement  fields  as 


H  ,  3n 
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(3.11) 


Using  (3-11)  and  the  definition  of  the  Stokes  corrections  it  is 
possible  to  write 
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The  residual  circulation  for  the  current  problem  is 
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and 
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(3.15) 


w*  =  w  +  i- 1-  YV 

N2  »y  ^ 


Using  the  eddy  thermodynamic  equation 


(a  +  uik)  -  f  v  '  ~  +  uikN2^  =  0 


and  multiplying  by  v'  and  averaging,  it  can  be  shown  that 

N2iku  v'c  -  fv'2  ~ 

TTY  - ^ - ~ -  (3.16) 

(a  +  iku) 

Using  (3-16)  and  comparing  (3.12)  and  (3.13)  with  (3.14)  and  (3.15) 
it  is  clear  that  the  Lagrangian-mean  velocity  and  the  residual 
velocity  are  equal  only  when  a  =  0.  When  a  =  0,  however,  both  the 
Lagrangian-mean  and  the  residual  circulation  are  zero  and  the 
transport  by  the  waves  and  the  mean  flow  exactly  counterbalance  each 
other.  Since  a  is  an  order  of  magnitude  or  more  smaller  than  ku 
below  the  sponge  layer  in  this  model,  the  Lagrangian-mean  and 
residual  circulations  should  be  similar.  When  transience  or 
external  heating  and  cooling  processes  are  allowed,  the  relationship 
between  the  two  circulations  becomes  more  complex.  In  Holton's 
(1981)  parameterization,  the  difference  between  the  Lagrangian-mean 
and  residual  circulations  is  contained  in  the  antisymmetric  part  of 
the  diffusion  tensor  and  is  assumed  to  be  negligible. 

In  the  next  section  the  Eulerian  velocities,  the  Lagrangian- 
mean  quantities,  and  the  residual  circulation  will  be  presented  for 
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steady,  dissipative  planetary  waves.  Zonal  mean  ozone  changes  can 
then  be  calculated  and  the  various  transport  methods  can  be 
compared . 

3.3  Results 

The  zonal  mean  ozone  distribution  used  in  the  steady  state 
model  is  shown  in  Fig.  3-1.  This  distribution  is  based  on  the 
October/November  values  given  by  Du'tsch  (1969)  and  is  similar  to  the 
distribution  used  by  Hartmann  and  Garcia  (1979). 

a.  Dynamic  quantities 

The  v  and  w  fields  are  shown  in  Fig.  3*2.  These  fields,  as 
expected,  show  rising  motion  in  the  northern  regions  of 
the  6-channel  and  sinking  motion  in  the  southern  regions.  The 
meridional  velocity  is  from  north  to  south.  This  flow  is  forced  by 
the  northward  eddy  heat  transport  of  the  planetary  waves. 

The  Stokes  drifts  as  calculated  with  (3.9b)  and  (3.9c)  are 
shown  in  Fig.  3.3.  These  drifts  can  be  viewed  as  the  effective 
velocity  fields  associated  with  the  wave.  Ihe  Stokes  drifts  are 
directed  in  an  opposite  sense  to  the  Eulerian  mean  velocity  and  are 
of  approximately  the  same  magnitude.  This  demonstrates  the  tendency 
for  the  wave  fields  to  counterbalance  the  mean  fields. 
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Figure  3.1.  Zonal  mean  ozone  distribution 
for  the  steady  state  model  (ppm). 
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Figure  3.2.  Eulerian  zonal  mean  velocities  (cm  sec-^). 
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The  Lagrangian-mean  velocity  fields,  which  are  the  sums  of  the 
Eulerian  zonal  mean  and  the  Stokes  drifts,  are  shown  in  Fig.  3.4. 
The  sense  of  the  circulation  is  rising  motion  in  the  southern 
regions  and  sinking  motion  in  the  polar  regions.  This  motion  is  the 
same  as  indicated  by  Schoeberl  (1981a).  Furthermore,  this  motion 
would  tend  to  enhance  the  diabatic  circulation  as  calculated  by 
Dunkerton  (1978).  The  Lagrangian-mean  velocity  is  of  much  smaller 
magnitude  than  either  the  Stokes  drifts  or  the  Eulerian  mean 
velocities . 

The  residual  circulation  is  shown  in  Fig.  3.5.  It  is 
geometrically  similar  to  the  Lagrangian-mean  velocity  but  the 
meridional  component  is  approximately  30%  larger  than  the 
Lagrangian-mean  meridional  velocity.  The  largest  differences 
between  the  two  circulations  occur  in  the  region  where  9u/3z  is 
largest.  Similar  increases  are  calculated  in  the  vertical  residual 
velocity  field. 
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Residual  velocities  (cm 


b.  Euler lan  transport 

The  continuity  equations  (3.7),  (3.8),  and  (3*10)  are  solved 
using  the  appropriate  dynamic  quantities.  Since  the  mean  chemistry 
appears  in  the  same  form  in  both  (3.8)  and  (3.10)  it  is  ignored  In 
order  that  it  might  not  obscure  the  differences  between  the  two 
formulations.  The  effects  of  the  mean  chemistry  term  will  be 
considered  separately. 

The  horizontal  (  vp  +  v^p"  )  and  vertical  (  wp  +  w'p>  ) 
fluxes  as  calculated  by  the  Eulerian  method  are  shown  In  Fig.  3.6. 
The  vertical  flux  Is  upward  in  the  northern  region  of  the  model  and 
downward  in  the  southern  region.  The  vertical  flux  is  confined 
primarily  to  the  conservative  region  and  reflects  the  Eulerian-mean 
vertical  velocity  field. 

The  horizontal  flux  field  is  quite  complicated  showing  two 
regions  of  northward  transport  separated  by  a  re*,  ion  of  southward 
transport .  Once  again  the  flux  is  largest  in  the  conservative 

region.  The  lower  region  of  northward  flux  is  dominated  by  eddy 
flux  as  described  by  Clark  and  Rogers  (1978).  The  region  of 

southward  flux  is  dominated  by  the  mean  meridional  velocity.  The 
upper  region  of  northward  transport  is  once  again  dominated  by  the 
eddy  flux.  In  the  upper  most  regions  of  the  model  both  the 

horizontal  and  vertical  fluxes  are  very  small  due  to  rapid 

photochemical  relaxation. 

The  rate  of  change  of  the  mean  ozone  field  is  calculated  as  the 
divergence  of  the  fluxes.  This  field  Is  shown  in  Fig.  3.7.  There 
is  an  increase  of  ozone  in  the  northern  section  and  a  decrease  in 
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the  southern  section.  Despite  the  fluxes'  being  largest  in  the 
conservative  region,  the  mean  field  change  is  concentrated  almost 
entirely  in  the  transition  region  in  agreement  with  the  results  of 
Hartmann  and  Garcia  (1979). 
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HORIZONTAL  FLUX  *  101  (ppm  m  sec  '» 


^  x  10~*  (ppm  ••c-1) 


Figure  3.7.  Mean  ozone  tendency,  ,  as  calculated 

O  t 

with  the  Eulerian  formulation  (ppm  sec-*). 
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c .  Lagrangian-mean  transport 

The  Lagrangian-mean  results  are  shown  in  Figures  3.8-3.10.  The 
advective  term  as  defined  in  Eq.  (3*10)  is  shown  in  Fig.  3.8.  In 
the  extreme  polar  region  there  is  a  decrease  of  ozone  above  the 
ozone  maximum  and  a  buildup  of  ozone  in  the  lower  stratosphere.  In 
the  southern  regions  there  is  a  decrease  in  the  lower  stratosphere 
and  an  increase  in  the  upper  stratosphere.  These  changes  are 
primarily  due  to  the  vertical  advection.  There  is  a  strong 
transport  of  ozone  northward  by  the  meridional  Lagrangian-mean 
velocity  in  the  midlatitudes.  The  Lagrangian-mean  velocity 
associated  with  planetary  waves  advects  ozone  poleward  and  downward. 

The  advective  term  calculated  using  the  residual  circulation 
shows  increases  and  decreases  in  the  mean  ozone  field  in  the  same 
regions  as  the  Lagrangian-mean  calculations.  In  all  regions  the 
magnitude  of  advection  by  the  residual  circulation  is  larger. 
Differences  as  large  as  30%  are  calculated  reflecting  the  increased 
magnitude  of  the  residual  circulation. 

Comparison  of  Figures  3-10  and  3.7  show  that  the  mean  field 
changes  are  the  same  when  calculated  by  the  Eulerian  or  Lagrangian- 
mean  methods.  Substituting  the  residual  velocity  for  the 
Lagrar.glan-mean  velocity  causes  no  significant  changes  in  the 
calculation  of  the  mean  field  changes . 


b.  residual  velocity 
Figure  3*8.  Mean  ozone  tendency 
caused  by  advection  (ppm  sec"*). 
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STIRRING  TERM  »  10  6  (ppm  sec  ') 


Y  (km) 

Figure  3-9-  Total  mean  ozone  tendency  as  calculated 
with  the  Lagrangian  formulation  (ppm  sec-1). 


d  ju 

x  10"®  (ppm  sec"') 


Figure  3.10.  Mean  ozone  tendency 
caused  by  "stirring”  (ppm  sec-1). 
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d.  The  "stirring"  term 


The  magnitude  of  the  stirring  term  is  shown  in  Fig.  3.9.  This 

field  quite  accurately  represents  the  ozone  time  rate  of  change 

shown  in  Fig.  3.10.  The  advection  by  the  Lagrangian-mean  velocity 
is  less  than  1%  of  the  stirring  term  in  regions  where  ozone  is 

photochemically  active.  The  stirring  term  is  identically  zero  when 
the  photochemistry  is  zero. 

In  order  to  understand  the  physical  meaning  of  the  stirring 
term  and  its  analogy  to  the  Eulerian  problem,  it  is  helpful  to 

compare  the  stirring  term  to  the  eddy  flux.  Assuming  that  the  wave 
is  steady  and  stationary  (3.7)  becomes 


1-iku  ^o 

A2  +  (ku)2  Ul  3xi 


(3.17) 


The  change  in  the  zonal  mean  ozone  field  due  to  the  divergence  of 
the  eddy  flux  (DEF)  is 
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The  stirring  term  (ST)  is 
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(3.19) 

ST  is  functionally  very  similar  to  DEF.  In  order  to  investigate  the 
precise  similarities  it  is  useful  to  consider  the  limits 
of  X  relative  to  the  dynamic  frequency  ku . 

For  X/ku  <  1, 


~  ~  =  — ~  Cl  "  -4-r)  (3-20) 

A*+(ku)2  (ku)2  (ku)2 


Using  (3-20),  (3-18)  and  (3-19)  can  be  rewritten  to  lowest  order  as; 


DEF  -  i  4r  Re<"PRe<uP  -S7 +  77:  ^ 

j  (ku)*-  J  i  (ku)  J  I 


(3.21) 


and 
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ST 


1  i_  P* 

P  3Xj  (ku)^ 


Re(u')Re(u')  ^ 


(3.22) 


Thus  for  small  A  the  divergence  of  the  eddy  flux  can  be  divided  into 
a  purely  chemical  and  purely  conservative  part .  The  conservative 
part,  the  second  term  on  the  RHS  of  (3.21)  is  the  conservative  flux 
described  by  Clark  and  Rogers  (1978).  The  chemical  part  of  DEF  is 
identical  to  ST;  therefore,  for  small  A  it  can  be  explicitly  shown 
that  the  stirring  term  represents  the  chemical  contribution  of  the 
divergence  of  the  eddy  flux. 

For  large  A  DEF  approaches  zero,  while  ST  approaches  some 
finite  value.  In  fact,  as  A  becomes  large  ST  loses  its  explicit 
chemistry  dependence  and  becomes 


ST  *  u 


s 

i 


(3.23) 


At  this  limit  both  the  Eulerian  mean,  (3.8),  and  Lagrangian-mean, 
(3.10)  continuity  equations  become 


at 


A  p 


(3.24) 


Eq.  (3.24)  indicates  that,  when  the  photochemistry  is  strong,  the 
wave  effects  become  small  so  that  the  mean  field  changes  are 
governed  by  the  mean  chemistry  and  the  advection  by  the  Eulerian- 
mean  velocity. 
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e .  Lagrangian  model 

While  the  analysis  above  reveals  the  relationship  of  the 
stirring  term  to  the  divergence  of  the  eddy  flux,  the  physical 
nature  of  the  stirring  term  remains  obscured  in  the  mathematics. 
Hartmann  (1981)  showed  how  pure  oscillatory  motion  of  a  fluid  parcel 
could  produce  changes  in  the  mean  ozone  density  in  a  fluid  tube.  A 
similar  Lagrangian  approach  can  be  used  here  to  investigate  the 
effects  of  the  planetary  wave  and  to  further  coalesce  the 
Lagranglan-mean  and  Eulerian  viewpoints. 

The  lagrangian  model  is  based  on  the  particle  trajectories  as 
determined  by  the  Lagrangian-mean  displacement  fields.  Using  (2.3) 
and  assuming  that  u  is  constant  and  F'  is  zero,  it  can  be  shown 
that  for  harmonic  disturbances  the  horizontal  displacement 
field,  n,  and  the  vertical  displacement  field,  C  are  given  by 

n  -  nQcos(az  +  kx)  (3.25) 

and 

C  -  Cqcos(o z  +  kx)  +  £isin(ctz  +  kx)  (3.26) 

where  a  is  the  vertical  wavenumber.  The  coefficients,  nQ»  CQ»  and 
are  functions  of  height,  meridional  distance,  and  the  dynamic 
parameters.  For  a  particular  value  of  y  and  z,  the  projection 
of  n  and  c  on  the  y,z  plane  is  an  ellipse  (Fig.  3.11).  These 
displacement  fields  are  similar  to  those  calculated  by  Matsuno 
(1980)  and  Danielsen  (1981). 
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The  displacement  fields  provide  a  model  of  the  parcel  motion 
associated  with  a  planetary  wave.  The  planetary  wave  displaces  a 
particular  parcel  from  its  equilibrium  position  and  causes  the 
parcel  to  move  in  a  helical  trajectory  around  the  globe.  The 
Lagrangian-mean  velocity  is  the  velocity  at  which  the  axis  of  the 
helix  moves.  For  the  small  amplitude,  steady  state  waves  of  this 
model  the  vertical  extent  of  the  ellipse  is  on  the  order  of  1-3  km 
and  the  horizontal  extent  is  on  the  order  of  hundreds  of  kilometers. 

First  consider  a  parcel  in  the  absence  of  chemistry.  The 
equilibrium  density  of  the  parcel,  pp,  is  uq(0,0).  If  the  parcel  is 
displaced  to  (n,C),  then  the  perturbation  density  p'  is  equal  to 
the  parcel  density  minus  the  background  density.  If  the 
displacements  are  small,  then 


Assuming  that  uq  is  a  function  of  z  only  and  that  it  increases  with 
height,  then  positive  ozone  perturbations  are  generally  associated 
with  northward  motion  and  negative  ozone  perturbations  are 
associated  with  southward  motion  (see  Fig.  3.11).  This  is 
equivalent  to  northward  eddy  flux  of  ozone.  Hence  this  is  the 
Lagrangian  description  of  the  conservative  flux  as  defined  by  Clark 
and  Rogers  (1978).  This  conservative  flux  is  largely  cancelled  by 
the  wave  induced  mean  circulation. 

In  the  case  when  5i  in  (3.26)  is  zero,  which  corresponds  to  an 
evanescent  wave,  the  trajectories  then  become  lines  instead  of 
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ellipses.  Ihen  the  perturbations  are  equally  correlated  with  both 


northward  and  southward  velocity  and  the  conservative  eddy  flux  is 
identically  zero. 


Figure  3.11.  Projection  of  displacement  fields 
n,  c  on  the  y,  z  plane  (km). 


Figure  3.12.  Parcel  orbiting  through  a  region  of 
discontinuous  photochemistry  4^  >  0. 

OZ 
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When  chemistry  Is  present  p'  still  changes  due  to  the  varying 
background  field,  but  the  chemistry  tries  to  bring  the  parcel's 
ozone  density  towards  the  background  value.  The  constituent  density 
for  the  parcel  is  determined  by 

-  S  -  Xp'  (3.27) 

where  S  represents  the  effective  source  caused  by  the  parcel  moving 
through  the  varying  background  density.  For  elliptical  trajectories 
S  is  a  harmonic  function  in  time,  proportional  to  the  magnitude  of 
the  velocity  and  the  background  constituent  gradients.  In  general, 
as  the  parcel  orbits  the  ellipse,  it  goes  through  varying 
photochemistry  so  that  X  is  also  a  function  of  time.  In  a  reference 
frame  moving  with  the  zonal  velocity  of  the  parcel,  (3.27)  can  be 
rewritten  as 


dp' 

dt"~ 


Sqcos  ut  -  X(t)p' 


(3.28) 


where  u  ■  ku. 

If  X  is  constant,  then  the  solution  to  (3.26)  for  large  t  can 
be  written  as 


S 

\t'  * - - -  (X  cos  u»t  +  u  sin  u»t)  (3.29) 

X2  +  u>2 


For  large  X,  p'  approaches  zero  which  is  equivalent  to  the  previous 
result  that  the  eddy  fluxes  disappear  in  the  fast  photochemical 
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region*  For  A  ■  0  the  solution  reduces  to  the  ozone  density 
appropriate  to  the  conservative  problem.  Note  that  the  average 
of  \i'  over  one  period  is  zero. 

In  the  regions  where  A  varies  strongly  over  the  extent  of  the 
ellipse,  the  perturbation  can  take  on  a  different  character.  Such 
regions  would  be  expected  to  exist  in  the  lower  parts  of  the 
transition  layer  and  at  the  polar  night  boundary  where  the  chemistry 
is  nearly  discontinuous.  Figure  3.12  shows  an  extreme  example  where 
the  chemistry  is  discontinuous  across  an  ellipse.  As  the  parcel 
travels  through  the  conservative  region  the  perturbation  density 
varies  as  S.  In  the  transition  region  (A/ku  ■  5)  the  perturbation 
is  decreased  by  ~  80%.  Assuming  the  background  density  Increases 
with  height ,  this  behavior  means  that  the  negative  perturbations  are 
selectively  damped  while  the  positive  perturbations  are  not. 
Therefore,  the  average  value  of  v'  over  one  period,  which  was  zero 
when  A  did  not  vary,  may  become  nonzero  in  the  presence  of  extreme 
chemical  gradients.  Such  regions  would  be  expected  to  have  a  great 
effect  upon  transport  calculations.  This  effect  is  analogous  to 
mean  density  changes  in  a  fluid  tube  as  calculated  by  Hartmann 
(1981). 
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Figure 


PARCEL  MIXING  RAV<0 


a.  constant  photochemistry 


TIME 


b.  photochemistry  at  the  bottom  of 
the  transition  layer 

3.13.  Solution,  at  large  t,  of  perturbation  ozone 
density  for  an  elliptically  orbiting  parcel. 
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Figure  3.13  demonstrates  this  effect  for  a  more  realistic 
problem.  Figure  3.13a  shows  (3.29),  the  steady  state  solution  to 
(3.28).  This  is  a  harmonic  function  with  zero  average.  Figure 
3.13b  shows  the  numerical  solution  to  (3.24)  for  an  ellipse  situated 
at  the  bottom  of  the  transition  layer  with  the  model 

photochemistry.  The  steady  state  solution  has  a  non-zero  time 
average,  as  expected. 

f.  The  effects  of  eddy  transport  on  the  mean  field 

Since  the  changes  in  the  mean  field  have  been  assumed  to  be 

O(a^),  the  calculation  of  the  steady  state  forcing  by  the  eddies  on 
the  mean  field  is  independent  of  the  mean  constituent  density. 

Therefore,  at  any  point  in  space  it  is  possible  to  write  the  mean 
field  equation  (3.8)  or  (3.10)  as 

=  -  Ap  +  F  (3.30) 

where  F  is  the  eddy  forcing  which  represents  the  divergence  of  both 
the  horizontal  and  vertical  fluxes.  F  is  independent  of  time; 

therefore,  the  solution  to  (3.30)  can  be  written  as 

u  -  x  (1  -  e"Xt)  (3-31) 

where  it  has  been  assumed  that  m»  the  change  in  the  equilibrium 
field,  is  zero  at  time  zero.  As  t  ♦  »  the  change  in  the  mean  field 
is  F/A. 
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In  the  conservative  region  the  mean  continuity  equation  can  be 
written  as  =  F  which  yields 


U  =  Ft  (3.32) 

In  order  to  investigate  the  effect  of  the  steady  state  waves  on  the 
mean  field  the  appropriate  equation  (3.31)  or  (3.32)  was  solved  at  t 
*  90  days.  This  calculation  would  represent  the  effects  of  the 
waves  if  they  persisted  for  the  entire  winter.  The  forcing  is  taken 
from  the  line  in  Fig.  3.7  through  the  region  of  maximum  ozone 
increase. 

Figure  3-14  shows  the  change  in  p.  The  large  rates  of  change 
that  are  calculated  in  the  transition  region  are  almost  completely 
counterbalanced  by  the  mean  chemistry.  At  the  altitude  of  maximum 
increase  in  Fig.  3.7  the  actual  change  in  the  mean  field  is  only  .5 
ppm  for  90  days. 

The  largest  changes  in  the  mean  field  are  in  the  lower 
stratosphere  where  changes  of  1.4  ppm  are  seen  at  90  days.  These 
changes  are  comparable  to  those  observed  in  the  northern 
hemisphere.  Immediately  below  the  transition  region  is  a  region  of 
depletion  caused  by  the  advection.  This  vertical  structure  of  the 
changes  in  the  mean  field  is  unrealistic  and  due  largely  to  the 
simplified  representation  of  photochemistry  and  dynamics  in  the 
lower  stratosphere. 
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HEIGHT  (km) 


CHANGE  IN  p  AT  90  DAYS 


Figure  3.14.  90  day  zonal  mean  ozone  change,  with  the 

mean  chemistry  included,  at  the  latitude 
of  the  maximum  ozone  increase  in  Fig  3.7  (ppm). 
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CHANGE  IN  p  AT  90  DAYS 


ppm 

a.  Southern 


CHANGE  IN  p  AT  90  DAYS 

60 

50 
E  40 

JC 

1  30 

uj 
I 

20 

10 
0 

4  -3-2-10  1  2 

ppm 

b.  Northern 

Figure  3.15.  90  day  zonal  mean  ozone  change,  with  the 

mean  chemistry  included,  at  the  northern 
and  southern  extremes,  as  determined  by 
Lagrangian-mean  advection  (LGM)  and 
residual  velocity  advections  (+)  (ppm). 
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The  relatively  large  positive  changes  in  the  mean  field  at  the 
bottom  of  the  transition  layer  are  in  the  region  where  the 


photochemistry  gradients  are  large.  This  is  the  region  where  a 
parcel  takes  on  a  non-zero  average  during  a  parcel  orbit.  The 
equilibrium  values  at  this  altitude  are  reached  after  5  days  of 
integration.  This  region  of  large  increase  just  above  a  region  of 
decrease  demonstrates  an  important  model  sensitivity  to  the  chemical 
structure . 

Since  the  advection  term  is  responsible  for  all  of  the  change 
in  the  mean  field  below  the  transition  region  it  is  worthy  of 
further  investigation.  Mean  field  tendencies  for  90  day  periods 
were  calculated  400  km  from  both  the  northern  and  southern 
boundaries  using  both  the  residual  and  tegrangian-mean  advection 
(Fig.  3.15).  The  large  changes  in  the  midlatitudes  caused  by 
meridional  transport  are  largely  negated  by  the  photochemistry. 

In  the  southern  region  the  change  in  the  mean  field  by 
advection  is  unlike  that  observed  in  the  stratosphere.  There  is  a 
depletion  of  nearly  4  ppm  at  24  km  and  an  enhancement  of  almost  2 
ppm  at  40  km.  In  the  northern  region  the  changes  in  the  lower 
stratosphere  are  not  unreasonable  showing  a  1.5  ppm  increase  at  24 
km.  However,  at  36  km  a  1  ppm  depletion  occurs.  It  is  obvious  that 
the  modeling  of  the  lower  stratosphere  is  inadequate.  Both  small 
scale  dynamics  and  photochemistry  which  may  act  to  counter  the 
tendencies  shown  here  have  been  neglected.  The  advection  calculated 
using  the  residual  circulation  produces  larger  changes  in  the  mean 
field  than  that  calculated  using  the  Lagrangian-mean  velocity.  This 
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means,  that  in  this  model,  the  antisymmetric  part  of  Holton's  (1981) 
diffusion  tensor  is  not  negligible. 

3.4  Conclusions 

The  changes  in  zonal  mean  ozone  as  forced  by  steady  state 
planetary  waves  have  been  calculated  using  both  a  Lagrangian-mean 
and  Eulerian  formulation.  Although  the  calculated  results  from  the 
two  methods  are  the  same,  the  use  of  the  two  methods  allows  for 
increased  interpretation  of  the  results. 

In  the  Eulerian  formulation  the  tendency  for  the  eddy  and  mean 
fluxes  to  cancel  is  well  illustrated.  Both  the  horizontal  and 
vertical  flux  fields  are  complex  and  largest  in  the  conservative 
region.  However,  the  changes  in  the  mean  field  in  the  absence  of 
the  mean  photochemistry  are  largest  in  the  transition  region  as 
found  by  Hartmann  and  Garcia  (1979). 

The  problem  of  compensation  by  two  large  terms  which  appears  in 
the  Eulerian  framework  also  appears  in  the  calculation  of  the 
Lagrangian-mean  velocity.  If  the  Lagrangian-mean  flow  could  be 
calculated  directly,  then  it  might  be  possible  to  avoid  the 
numerical  difficulties  of  compensation.  There  seems  to  be  no 
computational  advantage  of  the  Lagrangian-mean  method  to  calculate 
planetary  wave  transport  as  it  requires  the  formulation  of  Eulerian 
quantities,  their  conversion  to  Lagrangian-mean  quantities,  and  then 
the  calculation  of  the  transport. 

In  the  Lagrangian-mean  formulation  the  stirring  term  is 
dominant  in  the  absence  of  mean  chemistry.  The  stirring  is  equal  to 
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the  chemical  contribution  of  the  divergence  of  the  eddy  flux 
when  \/uk  <  1.  When  X  is  large  the  stirring  term  is  equivalent  to 
the  Stokes  advection.  Using  particle  displacement  fields  it  is 
possible  to  predict  the  type  of  eddy  transport  that  is  to  be 
expected  in  a  certain  region.  It  was  found  that  \i'  has  a  non-zero 
time  average  in  the  presence  of  strong  chemical  gradients. 
Calculation  of  the  mean  field  changes  shows  that  the  mean  chemistry 
is  able  to  compensate  for  the  eddy  chemistry  changes  except  where 
the  chemical  gradients  are  large.  This  indicates  that  the 
interaction  of  steady  planetary  waves  and  photochemistry  need  only 
be  modeled  in  the  presence  of  large  chemical  gradients. 

For  seasonal  integrations  the  largest  changes  in  the  mean  field 
are  caused  by  advection  by  the  Lagrangian-mean  velocity  in  the 
conservative  region.  The  advection  is  so  weak  that  in  the  chemical 
regions  its  effects  are  negligible.  Using  the  residual  circulation 
instead  of  the  Lagrangian-mean  circulation  overestimates  the 
advection  by  the  planetary  waves.  In  this  particular  model  the 
largest  differences  occur  between  the  Lagrangian-mean  and  residual 
circulations  in  chemically  active  regions.  Hence,  the  difference 
between  the  ozone  field  changes  calculated  with  Lagrangian-mean  and 
residual  circulations  is  small.  This  would  not  necessarily  be  true 
for  constituents  with  different  chemistry. 

The  transport  mechanisms  investigated  here  are  probably  not 
adequate  to  model  planetary  wave  effects  in  the  stratosphere. 
Sudden  warmings,  which  are  highly  transient,  will  certainly  have  a 
great  effect  upon  the  ozone  distribution  (Garcia  and  Hartmann, 
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1980).  Therefore,  steady  state  planetary  wa .  e  transport  models  are 
probably  only  relevant  in  the  late  fall  and  early  winter.  The 
planetary  wave  transport  would  have  to  be  supplemented  by  transport 
due  to  the  diabatic  circulation  and  small  scale  dynamics  in  the 
lower  stratosphere.  Transport  by  planetary  waves  during  sudden 
warmings,  transport  by  the  diabatic  circulation,  and  the  relative 
importance  of  the  two  circulations  will  be  presented  in  the  next 
chapter. 
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Chapter  4 


TRANSPORT  BY  TRANSIENT  PLANETARY  WAVES  AND  THE  DIABATIC  CIRCULATION 

4.1  Introduction 

In  this  chapter  the  transport  of  ozone  by  time  dependent 

planetary  waves  on  a  8-plane  is  discussed.  In  Section  4.2  the 

results  of  the  dynamic  model  for  both  600  m  and  900  m  topographic 
forcing  are  presented.  In  Sections  4.3,  4.4,  and  4.5  the  transport 
of  ozone  by  the  diabatic  circulation  alone,  by  planetary  waves 

alone,  and  by  the  combination  of  the  diabatic  and  planetary  wave 
circulations  is  discussed .  The  results  are  analyzed  in  Section  4.6. 

4.2  Dynamics 

Equations  2.3  and  2.4  are  solved  on  a  3-plane  centered  at  60°  N 
for  zonal  wavenumber  one  perturbations.  The  wave  is  forced  by  zonal 

flow  over  a  wavenumber  one  topography  with  a  height  of  either  600  m 

or  900  m.  With  forcing  of  this  magnitude  there  is  at  least  one 

sudden  warming  accompanied  by  a  reversal  in  the  zonal  mean  wind 

field.  For  long  Integrations  a  steady  state  solution  is  reached 
with  permanent  easterlies  at  some  level.  Schoeberl  (1982b)  presents 
a  detailed  analysis  of  the  model  dynamics. 

During  winter  it  is  common  to  see  distinct  pulses  of  wave 
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amplitude  and  either  major  or  minor  warmings  (Labltzke  and  Goretzki, 
1982).  After  a  major  or  minor  warming  the  stratosphere  usually 
returns  to  its  winter  climatological  state.  If  the  warming  occurs 
late  in  the  winter,  then  the  zonal  mean  circulation  after  the 
warming  evolves  into  the  summertime  regime  where  easterlies  are 
present  throughout  the  entire  stratosphere.  This  type  of  warming  is 
called  a  final  warming.  The  dynamical  model  used  here  does  not 
attempt  to  represent  all  of  the  features  of  the  winter  stratosphere 
only  to  simulate  warming  events  and  their  corresponding  circulation 
patterns . 

The  details  of  the  numerical  scheme  were  given  in  Section 
2.3b.  Results  will  be  presented  as  time-height  cross  sections  at 
the  center  of  the  channel  and  at  a  northern  and  southern  point, 
corresponding  to  60°  N,  46°  N,  and  74°  N  on  the  8-plane. 

a.  600  m  topography 

Time  series  of  60  days  for  u,  the  zonal  mean  wind,  and  the 
magnitude  of  $  ,  the  planetary  wave  geopotential,  are  shown  at  the 
channel  center  in  Figures  4.1  and  4.2  for  the  600  m  orography.  The 
wave  amplitude  increases  at  all  altitudes  until  day  22  when  it 
reaches  a  maximum  of  about  1.6  gpkm.  After  day  22  the  wave 
amplitude  decreases  rapidly  until  day  30  and  continues  to  decrease 
more  slowly  until  day  47.  On  day  47  the  wave  amplitude  starts  to 
increase  at  high  altitudes  and  is  still  Increasing  at  the  end  of  the 
time  series. 
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;ure  4.1.  Zonal  mean  wind  for  600  m  topography 
at  the  channel  center  (m  sec-1). 
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Figure  4.2.  Geopotential  amplitude  for  600  m  topography 
at  the  channel  center  (gpdm). 
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Coincident  with  the  maximium  in  the  wave  amplitude  field  is  the 
appearance  of  easterlies  in  the  zonal  mean  wind  field.  The  presence 
of  easterlies  prohibits  further  vertical  propagation  of  the 
stationary  wave  (Charney  and  Drazin,  1961)  and,  hence,  the  rapid 
decrease  of  the  wave  amplitude.  The  easterly  zonal  wind  reaches  a 
maximum  velocity  of  53  m  sec-^  on  day  29  and  then  decreases  in 
intensity  until  day  45  when  westerlies  begin  to  appear  again  at  the 
edges  of  the  channel.  The  reappearance  of  westerlies  permits 
vertical  planetary  wave  propagation  subsequently  leading  to  the 
second  increase  in  wave  amplitude  at  day  47. 

Since  it  has  been  predicted  that  transport  in  the  vicinity  of 
critical  lines  might  be  large  (Schoeberl,  1981b;  Matsuno  and 
Nakamura,  1979),  the  descent  of  the  zero  wind  line  (critical  line 
for  stationary  waves)  is  of  particular  importance  to  the  study  of 
tracer  transport.  For  the  600  m  topography  the  critical  line  in  the 
center  of  the  channel  has  descended  to  48  km  on  day  25  and  reaches  a 
minimum  altitude  of  36  km  on  day  35.  The  critical  line  remains  in 
the  neighborhood  of  36  km  and  then  starts  to  rise  on  day  40.  This 
rise  is  related  to  the  reappearance  of  westerlies  at  high  and  low 
latitudes.  The  critical  line  ascends  to  50  km  on  day  45  and  then 
starts  to  descend  again  as  the  wave  amplitude  increases. 
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ure  4.3.  Zonal  mean  wind  for  900  m  topography 
at  the  channel  center  (m  sec-*). 
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Figure  4.4.  Geopotential  amplitude  for  900  m  topography 
at  the  channel  center  (gpdm). 
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b.  900  o  topography 

The  mean  zonal  wind  field  and  the  wave  amplitude  field 
associated  with  the  900  m  topography  are  shown  In  Figures  4.3  and 
4.4  for  a  120  day  Integration.  These  fields  are  similar  in  many 
ways  to  the  600  m  topography  fields.  The  wave  amplitude  Increases 
rapidly  until  about  day  20  where  it  reaches  a  maximum  of  1.6  gpkm. 
The  fact  that  the  forcing  has  been  Increased  by  50%  but  the  wave 
amplitude  reaches  the  same  maximum  amplitude  is  consistent  with  the 
wave  saturation  theory  of  Schoeberl  (1982a).  The  main  effect  of  the 
increased  forcing  is  to  cause  the  wave  maximum  to  be  reached 
somewhat  earlier.  The  wave  amplitude  rapidly  decreases  between  days 
20  and  30  after  which  it  slowly  decreases  until  day  50.  There  is 
another  rapid  increase  in  wave  amplitude  after  day  50  which 
maximizes  on  day  65.  The  wave  amplitude  decrease  following  this 
maximum  is  not  as  rapid  as  that  following  the  initial  warming  event. 

Coincident  with  the  first  wave  maximum  is  the  rapid 
deceleration  of  the  mean  zonal  wind  and  the  formation  of  easterlies 
which  reach  a  maximum  velocity  on  day  26  of  57  m  sec-^.  The  second 
wave  maximum  on  day  65  is  not  instantaneously  associated  with  the 
onset  of  high  level  easterlies.  The  wave  forcing  of  the  mean  flow 
from  days  55-65  decelerates  the  mean  zonal  wind  restoring  the  low 
level  easterlies  which  then  prohibit  propagation  of  the  wave  into 
the  stratosphere.  The  low  level  easterlies  ultimately  disappear 
about  day  100  which  allows  the  wave  to  once  again  propagate 
vertically.  The  zonal  mean  westerlies  are  so  weak  at  this  time  that 
even  the  relatively  weak  forcing  associated  with  the  wave  pulse  at 


83 


day  100  Is  capable  of  producing  easterlies  throughout  a  deep  portion 
of  the  atmosphere. 


The  zero  wind  line  reaches  minimum  altitudes  of  26  km  on  day  37 
and  24  km  on  day  76.  Between  days  45  and  55  the  critical  line  rises 
and  the  wave  propagates  and  rebuilds  the  low  level  easterlies.  The 
low  level  critical  line  is  absent  for  a  brief  period  beginning  at 
day  100 . 

c.  The  diabatic  circulation 

As  mentioned  in  the  introduction,  early  theories  of 
stratospheric  circulation  postulated  a  direct  cell  with  rising 
motion  in  the  equatorial  region  and  sinking  motion  in  the  polar 
region.  This  circulation  was  thought  to  describe  the  dryness  of  the 
stratosphere  and  the  buildup  of  ozone  at  high  latitudes  (Brewer, 
1949;  Dobson,  1956).  Observations  during  the  1960 's  showed, 
however,  that  the  northern  hemisphere  stratospheric  mean  circulation 
consisted  of  two  cell6  with  rising  motion  over  the  pole  and  the 
equator.  Dunkerton  (1978)  pointed  out  that  in  the  absence  of  waves 
the  diabatically  driven  zonal  mean  field  is  a  good  approximation  to 
the  Lagrangian-mean  velocity  field  and  is  consistent  with  the 
Brewer-Dobson  circulation.  Hence,  the  large  scale  mass  transport 
associated  with  radiative  forcing  is  described  by  a  one  cell 
circulation  as  originally  postulated. 

The  diabatic  circulation  used  here  is  derived  from  Dunkerton 
(1978)  and  the  more  detailed  calculations  of  Murgatroyd  and 
Singleton  (1961).  The  circulation  appropriate  for  the  winter 


solstice  is  shown  in  Fig.  4.5.  The  time  variation  is  assumed  to  be 
harmonic  such  that  the  intensity  of  the  flow  reaches  a  maximum 
shortly  after  the  winter  solstice  and  then  decreases  and  reaches  a 
minimum  in  the  spring. 

d.  Initialization 

The  zonal  mean  ozone  field  used  to  initialize  the  transport 
model  is  shown  in  Fig.  4.6.  This  field  is  derived  from  the  data  of 
Diitsch  (1969)  and  is  representative  of  mid-November.  Data  north  of 
70°  and  above  45  km  are  scarce  and  of  low  quality.  In  these  regions 
the  prescribed  ozone  field  is  based  on  the  vertical  profile  derived 
by  Kreuger  and  Minzner  (19  76)  and  the  ozone  field  used  by  Fels  ej^ 
al .  (1980).  The  eddy  ozone  field  is  initially  assumed  to  be  zero. 

A  review  of  stratospheric  warming  data  indicates  that  there  is 
either  a  major  or  minor  warming  during  the  first  or  second  week  of 
January.  Since  the  model  warmings  occur  about  three  weeks  after  the 
planetary  waves  are  turned  on,  December  15  has  been  chosen  as  the 
initialization  dat_.  In  order  to  properly  correlate  the  ozone  field 
with  the  planetary  wave  dynamics  the  Initial  ozone  field  is  advected 
with  the  diabatic  circulation  until  Dec  15. 
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Figure  4.5.  Diabetic  circulation 
at  the  winter  solstice  (cm  sec“^). 
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4.3  Transport  by  the  diabatic  circulation 

The  change  in  the  mean  ozone  field  (anomaly)  due  solely  to  the 
diabatic  circulation  is  presented  in  one  year  time  series  in  Figures 
4.7  and  4.8.  In  general  there  is  a  positive  anomaly  in  the  lower 
stratosphere  (20-30  km)  for  much  of  the  year  and  a  negative  anomaly 
in  the  autumn.  Since  the  ozone  is  in  near  equilibrium  above  35  km, 
a  positive  (negative)  anomaly  is  associated  with  an  increase 
(decrease)  in  total  ozone.  Using  this  correlation  it  is  seen  that 
the  total  ozone  reaches  a  maximum  in  early  April  and  a  minimum  in 
October  as  opposed  to  an  observed  March  maximum  and  October  minimum. 

There  is  a  strong  latitudinal  dependence  of  the  anomaly.  The 
wintertime  increase  is  larger  and  reaches  lower  altitudes  at  high 
latitudes.  In  the  northern  region  (Fig.  4.8)  a  maximum  increase  of 
3.5  ppm  is  centered  around  26  km,  and  increases  greater  than  1  ppm 
are  seen  as  low  as  15  km.  At  the  center  of  the  channel  (Fig.  4.7) 
an  anomaly  of  3  ppm  is  centered  at  30  km,  wMle  in  the  southern 
region  (not  shown)  a  maximum  increase  of  .5  ppm  is  centered  around 
32  km.  The  largest  increase  in  the  entire  domain,  4.25  ppm,  occurs 
at  70°. 

The  computed  total  ozone  field  (Fig.  4.9)  differs  considerably 
from  observations.  The  time  series  for  Dec . 15-Apr . 15  shows  a 
maximum  of  538  DU,  an  overestimate  of  20%,  which  is  reached  about 
the  first  of  April  at  70°.  Observations  of  ozone  in  the 
stratosphere  (Fig.  2.11)  show  a  maximum  near  the  pole  about  March  1 
of  approximately  450  DU. 
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Figure  4.9.  Latitude-time  cross  section  of  total  ozone 
as  calculated  with  the  diabetic  circulation  (DU). 
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The  changes  in  ozone  above  45  km  in  Figures  4.7  and  4.8  are 
negligible  due  to  the  strong  photochemistry.  Therefore,  where  the 
diabatic  circulation  is  strongest  no  changes  in  the  ozone  field  are 
found.  The  transport  effects  by  the  diabatic  circulation  are  most 
notable  in  the  transition  and  conservative  regions  and  are  due 
primarily  to  the  vertical  velocity  which  is  downward  at  all 
latitudes  and  constantly  advecting  ozone  into  the  lower 
stratosphere.  The  maximum  anomaly  occurs  about  500  km  north  of  the 
maximum  vertical  velocity  indicating  the  relative  ineffectiveness  of 
the  meridional  velocity  (<.3  m  sec-*)  at  transporting  ozone 
northward . 

Because  insolation  is  weak  during  the  winter  the  ozone 
photochemistry  is  much  slower  than  modeled  here  at  high  latitudes. 
The  overestimation  of  the  photochemical  source  term  may  be 
responsible  for  the  unrealistic  buildup  in  the  northern  region.  The 
results  from  two  experiments  where  latitudinal  variation  of  the 
chemistry  is  included  are  shown  in  Fig.  4.10.  In  the  first 
experiment  the  photochemistry  is  given  a  cosine  of  latitude 
dependence  such  that  the  chemical  relaxation  weakens  and  disappears 
at  the  northern  boundary.  This  is  the  same  latitudinal  dependence 
as  used  by  Garcia  and  Hartmann  (1980).  Comparing  Figures  4.8  and 
4.10a  indicate  there  is  no  significant  difference  in  the  results 
from  the  case  with  no  latitudinal  variation  of  the  chemistry. 
Apparently,  the  diabatic  circulation  is  so  weak  that  even  with  much 
slower  photochemistry  the  ozone  never  varies  much  from  its 


equilibrium  value . 


In  the  second  experiment  (Fig  4.10b)  the  photochemistry  is  set 
exactly  equal  to  zero  north  of  67.5°  for  the  entire  time 
Integration.  The  positive  anomaly  in  the  lower  stratosphere  remains 
the  same  as  in  the  other  experiments,  but  there  is  a  negative 
anomaly  around  50  km  of  more  than  1  ppm.  This  region  of  depletion 
is  due  to  the  Inability  of  the  ozone  to  return  to  its  equilibrium 
state.  Since  the  Increase  in  the  lower  stratosphere  remains  the 
same,  some  other  process  besides  the  latitudinal  variation  of  the 
photochemistry  must  be  responsible  for  prohibiting  the  excessive 
buildup  of  ozone  at  70° . 

In  summary,  the  dlabatic  circulation  during  winter  is  directed 
such  that  transport  is  poleward  and  downward.  The  northward 
component  of  the  velocity  is  too  small  to  effectively  transport 
ozone  to  polar  regions;  however,  the  vertical  component  acts  over  a 
considerable  depth  of  the  stratosphere  and  transports  significant 
amounts  of  ozone  into  the  lower  stratosphere.  The  total  ozone 
calculated  using  the  dlabatic  circulation  alone  is  too  large  in 
magnitude.  In  addition,  the  maximum  amounts  occur  too  fai  south  and 
about  1  month  too  late . 
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Figure  A. 10.  Change  in  ozone  at  74°  with  horizontally 
varying  chemistry  (ppm). 
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4.4  Planetary  wave  transport 

Transport  by  steady,  stationary,  dissipative  planetary  waves 
was  discussed  in  Chapter  3.  It  was  found  that  in  the  conservative 
region  of  the  stratosphere  planetary  waves  tended  to  advect  ozone 
from  low  latitudes  towards  the  pole.  In  the  transition  region  it 
was  found  that  stirring  of  parcels  by  the  waves  across  regions  of 
widely  varying  photochemistry  could  produce  changes  in  the  mean 
field  that  the  mean  chemistry  could  not  completely  counterbalance. 
In  this  section  the  results  of  transport  by  transient  planetary 
waves  will  oe  presented.  In  order  to  study  transport  mechanisms  it 
is  useful  to  study  both  conservative  transport  (chemically  inactive) 
and  chemical  transport. 

a.  Conservative  transport 

The  changes  in  the  mean  ozone  field  for  conservative  transport 
by  planetary  waves  forced  with  600  »  orography  is  shown  in  Figures 
4.11  and  4.12  for  the  southern  and  northern  points.  Centered  around 
40  km  is  a  zone  of  poleward  ozone  transport  and  just  above  this 
zone,  centered  at  50  km,  is  a  region  of  equatorward  transport.  With 
the  900  m  orography  (Fig.  4.13),  similar  structure  is  evident  with 
the  northward  transport  centered  at  35  km  and  the  southward 
transport  centered  at  50  km.  The  maximum  positive  change  at  the 
northern  point  for  the  900  m  case  is  2.2  ppm  and  occurs  on  day  32; 
for  the  600  m  case  a  maximum  of  slighlty  more  than  1  ppm  occurs  on 
day  35.  For  both  cases  there  is  relatively  little  change  at  the 
center  of  the  channel. 
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Figure  4.11.  Change  In  ozone  for  600  m  topography 
at  46°  with  no  chemistry  (ppm). 
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Figure  4.12.  Same  as  Fig.  4.11  but  at  74°. 


Figure  4.13.  Same  as  Fig.  4.12  but  for  900  m  topography 
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Superimposed  on  Figures  4. 11-4. 13  Is  the  location  of  the 
critical  line  at  the  appropriate  latitude.  For  all  cases  the  level 
of  maximum  northward  transport  is  coincident  with  the  minimum 
descent  of  the  critical  line.  The  layer  of  southward  transport  is 
above  the  critical  line.  The  transport  calculated  here  supports  the 
theory  of  Schoeberl  (1981b)  which  predicts  that  a  wave  energy 
absorbing  critical  line  will  produce  a  jet  of  northward  transport  at 
the  critical  line.  This  transport  mechanism  will  be  discussed 
further  in  Section  4.6. 

b.  Chemical  transport 

The  transport  of  ozone  for  both  600  m  and  900  m  orography  with 
photochemistry  is  presented  in  Figures  4.14-4.16.  There  is  the 
general  tendency  for  the  upper  region  of  equatorward  transport  to  be 
significantly  reduced  in  all  cases  due  to  the  photochemistry.  With 
photochemistry  the  transport  appears  to  occur  in  a  thin  layer  in  the 
conservative  region  and  at  the  bottom  of  the  transition  region. 

For  the  600  m  orography  at  the  southern  point  (Fig  4.14)  the 
magnitude  of  the  depletion  at  40  km  has  decreased  by  35%  with  the 
addition  of  photochemistry  (compare  to  Fig.  4.11).  This  might  be 
expected  since  the  chemistry  would  tend  to  moderate  large 
differences  from  equilibrium.  The  magnitude  of  the  positive  anomaly 
at  the  northern  grid  point  (Fig.  4.15)  on  the  other  hand  has  not 
moderated  with  the  addition  of  chemistry  but  has  increased  by  14%. 

At  the  northern  point,  with  the  900  m  orography  (Fig  4.16),  the 
anomaly  field  is  once  again  more  jet  like  than  without  chemistry. 
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but  the  magnitude  has  remained  the  same.  The  transport  calculated 
at  the  southern  point  shows  similar  behavior  with  the  low  level 
negative  anomaly  maintaining  the  same  magnitude  with  the  addition  of 
photochemistry.  The  level  of  maximum  northward  transport  for  the 
900  m  case  occurs  significantly  lower  in  the  stratosphere  than  with 
the  600  m  forcing;  therefore,  the  addition  of  chemistry  would  be 
expected  to  have  relatively  little  effect  on  the  transport. 

It  is  of  interest  to  examine  the  extended  time  series  of  the 
900  m  forcing  (Fig.  4.16).  In  this  figure  the  solid  heavy  line  is 
the  critical  line  at  the  latitude  of  the  figure  and  the  heavy  dashed 
line  is  the  critical  line  at  the  center  of  the  channel.  While  the 
critical  line  at  high  and  low  latitudes  ascends  and  dissipates, 
there  is  still  intense  south  to  north  transport  centered  around  35 
km.  This  transport  is  obviously  associated  with  the  critical  line 
at  the  center  of  the  channel.  The  anomaly  decreases  and  increases 
coincident  with  the  ascent  and  descent  of  the  critical  line. 

The  largest  northern  anomaly  during  the  entire  time  integration 
occurs  on  day  70  when  the  critical  line  reaches  its  lowest  altitude 
at  the  center  of  the  channel.  At  this  time  there  is  no  critical 
line  present  at  any  altitude  at  the  northern  point.  This  result 
shows  that  the  horizontal  extent  of  the  transport  associated  with 
the  critical  line  at  the  center  of  the  channel  is  significant. 
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Figure  4.16.  Same  as  Fig.  4.15  but  for  900  m  topography, 
(bold  dashed  line  is  the  critical  line  at  the  center  of 
the  channel) 
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4.5  Transport  by  both  waves  and  the  diabatic  circulation. 

The  change  in  the  zonal  mean  ozone  field  associated  with  the 
diabatic  circulation  superimposed  on  the  900  m  topography  dynamics 
is  shown  in  Figures  4.17-4.19.  The  900  m  topography  was  chosen  for 
these  experiments  since  this  amplitude  forcing  is  more  consistent 
with  the  terrestrial  topography  (Berkofsky  and  Bertoni,  1955),  and 
the  critical  line  for  this  case  descends  below  32  km.  One  of  the 
WMO  criteria  for  the  classification  of  a  major  stratospheric  warming 
is  that  the  wind  reversal  descends  below  10  mb  which  is 
approximately  32  km.  The  transport  by  the  combined  circulation  is 
much  different  from  the  transport  by  each  individual  type  of 
circulation  indicating  that  both  the  wave  and  diabatic  transport 
regimes  are  important  in  the  lower  stratosphere. 

At  the  center  of  the  channel  (Fig.  4.18)  there  is  a  maximum 
increase  between  25  and  30  km  of  1.25  ppm  in  late  March.  There  is  a 
series  of  short  term  decreases  around  40  km  with  the  most  intense 
decrease  of  1.25  ppm  also  occurring  in  late  March.  The  region  of 
steady  increase  is  due  to  the  downward  transport  by  the  diabatic 
circulation.  The  intermittent  nature  of  the  ozone  depletion  at  40 
km  is  associated  with  the  wave  transport.  With  the  combined 
circulation  the  increase  in  ozone  at  the  channel  center  is  reduced 
by  more  than  50%  when  compared  to  the  changes  caused  by  the  diabatic 
circulation  alone  (Fig.  4.7),  which  indicates  that  the  planetary 
waves  are  transporting  ozone  out  of  the  center  of  the  channel. 
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At  the  northern  latitude  (Fig.  4.29)  there  is  a  uniform 

increase  of  ozone  at  all  levels  with  the  maximum  increase  of  3.12 
ppm  occurring  close  to  30  km,  about  4  km  lower  than  the  maximum 
associated  with  the  planetary  waves  alone  (Fig.  4.16)  and  4  km 
higher  than  that  associated  with  the  diabatic  circulation  alone 
(Fig.  4.8).  An  anomaly  of  1  ppm  is  evident  at  25  km,  which  is  10  km 
higher  than  a  similar  anomaly  calculated  using  on]y  the  diabatic 
circulation  (compare  to  Fig.  4.8).  The  largest  anomaly  occurs  in 
early  March  and  is  coincident  with  the  lowest  altitude  of  the 
critical  line.  Following  the  maximum  in  early  March  the  positive 
anomaly  remains  relatively  stable  at  a  maximum  value  somewhat 
greater  than  2  ppm. 

At  the  southern  latitude  (Fig.  4.17)  there  is  a  region  of 
relatively  strong  poleward  transport  between  35  and  40  km  and  below 
this  layer  is  a  region  of  increase.  As  at  the  center  of  the 

channel,  the  region  of  increase  is  associated  with  the  diabatic 
circulation  and  the  decrease  is  associated  with  the  planetary 

wave.  Such  large  depletions  are  not  observed  at  this  latitude  (44°) 
in  the  stratosphere. 

The  total  ozone  time  series  (Fig.  4.20)  shows  much  better 
agreement  with  observation  at  high  latitudes  than  the  total  ozone 
field  calculated  with  the  diabatic  circulation  (Fig  4.9).  A  late 

winter  maximum  of  almost  440  DU  occurs  during  early  March  at  the 
most  northern  latitudes.  There  is  an  earlier  maximum  of  slightly 
more  than  420  DU  coincident  with  the  initial  descent  of  the  critical 
line.  This  earlier  maximum  produces  a  wave  like  feature  that  is 
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similar  to  the  wave  structure  shown  in  Fig  2.10.  Both  maxima  are 
associated  with  the  critical  line  descent  after  the  first  warming 
event  (Fig.  4-3). 

In  the  southernmost  section  of  the  model  there  is  too  little 
ozone  from  the  beginning  of  the  combined  planetary  wave-diabatic 
Integration  (Dec.  15).  Furthermore,  during  the  initial  descent  of 
the  critical  line  there  is  excessive  transport  out  of  the  south. 
The  excessive  transport  occurs  at  the  time  that  the  critical  line 
descends  to  low  altitudes  at  the  northern  and  southern  points. 


105 


4.6  Discussion. 


The  transports  associated  with  the  diabatic  and  planetary  wave 
circulations  are  fundamentally  different  from  each  other  and  of 
equal  importance.  The  diabatic  circulation  is  primarily  responsible 
for  transport  of  ozone  downward  from  the  photochemical  source  region 
while  the  planetary  wave  circulation  is  responsible  for  north-south 
transport.  The  combined  circulations  do  a  particularly  good  job  of 
reproducing  the  features  of  the  total  ozone  distribution  at  high 
latitudes . 

The  vertical  structure  of  the  winter  ozone  field  is  not  as  well 
simulated  as  the  total  ozone  distribution.  The  computed  increase  of 
more  than  3  ppm  at  high  latitudes  is  concentrated  at  30  km  which  is 
larger  than  observed  at  this  altitude.  The  increase  between  20  and 
25  km  is  significantly  smaller  than  is  observed.  The  increase  in 
ozone  at  high  latitudes  is  due  largely  to  the  downward  transport  of 
ozone  out  of  the  photochemical  region  by  the  diabatic  circulation  at 
70°  and  the  subsequent  northward  transport  of  this  ozone  to  high 
latitudes  by  the  planetary  waves. 

Enough  ozone  is  transported  out  of  the  source  region  to  high 
latitudes  to  model  the  total  ozone  maximum;  however,  redistribution 
of  ozone  within  the  conservative  region  is  needed  to  accurately 
reproduce  the  ozone  distribution.  Synoptic  scale  disturbances  are 
probably  instrumental  in  determining  the  vertical  structure  of  ozone 
within  the  lower  stratosphere. 

In  the  southern  region  of  the  channel  and  to  a  lesser  degree  in 
the  middle  of  the  channel  the  vertical  structure  is  more  poorly 
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represented.  There  is  excessive  transport  out  of  the  south  by  the 
planetary  waves  between  35  and  40  km  which  leads  to  unrealistic 
vertical  profiles.  This  excessive  transport  is  also  evidenced  by 
the  underestimation  of  total  ozone  at  low  latitudes  during  mid- 
January.  The  overestimation  of  transport  out  of  the  southern  region 
might  be  due  to  an  Inadequate  representation  of  chemistry  at  low 
altitudes;  the  underestimation  of  the  transport  by  the  diabatic 
circulation  into  the  south  from  either  high  altitudes  or  lower 
latitudes;  geometric  effects  of  the  0-plane  transport;  and/or  to  the 
neglect  of  small  scale  dynamics. 

In  order  to  understand  the  particular  nature  of  the  transport, 
it  is  necessary  to  further  analyze  the  transport  mechanisms 
involved.  It  was  mentioned  in  Chapter  3  that  even  for  steady, 
stationary,  dissipative  planetary  waves,  there  were  computational 
difficulties  in  calculating  lagrangian-mean  quantities.  In  the 
presence  of  highly  transient  dynamics  these  difficulties  are  more 
severe,  and  the  calculation  of  Lagrangian-mean  velocities  is  time 
consuming  with  the  result  being  uncertain  (McIntyre,  1980). 
Furthermore,  in  the  vicinity  of  the  critical  lines,  the  displacement 
vectors  are  nearly  singular  and  their  calculation  is  difficult. 
Ihmkerton  et  al .  (1981)  also  experienced  difficulty  in  the 

calculation  of  Lagrangian-mean  quantities  with  realistic  dynamical 
models.  Therefore  the  transport  will  be  analyzed  in  the  Eulerian 
framework  with  comments  made  on  the  Lagrangian-mean  interpretation 
where  it  is  appropriate. 
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a.  Diabatic  circulation 

At  all  times  during  the  winter  at  extratropical  latitudes  the 
flux  due  to  the  diabatic  circulation  is  northward  and  downward.  The 
northward  transport  is  quite  weak  and  does  little  to  change  the 
ozone  distribution.  The  vertical  transport  is  effective  at 
transporting  ozone  from  the  photochemical  and  transition  regions 
into  the  upper  part  of  the  conservative  region. 

While  the  diabatic  circulation  does  transport  ozone  northward 
and  downward  it  is  inadequate  alone  to  explain  the  wintertime  polar 
ozone  buildup.  The  meridional  component  of  the  diabatic  circulation 
is  large  enough  to  cause  significant  transport  only  at  high 
altitudes  where  the  photochemistry  is  dominate.  Therefore  Brewer- 
Dobson  circulations,  which  involve  the  upward  transport  of  ozone  at 
low  latitudes;  northward  transport  at  high  altitudes  in  the 

midlatitudes;  and  subsidence  at  high  latitudes,  cannot  explain  ozone 
transport . 

b.  Planetary  waves 

The  planetary  wave  transport  is  more  complex  than  that 

presented  in  Chapter  3  for  steady-state  transport.  Examination  of 
Fig.  4.16,  for  Instance,  shows  that  the  transport  is  highly  time 

dependent  with  strong  pulses  of  northward  transport.  The  planetary 
wave  transport  was  demonstrated  to  be  directly  related  to  the 

descent  of  the  critical  line  in  the  previous  section. 

According  to  theory  (Charney  and  Drazin,  1961),  it  is 

impossible  for  stationary  planetary  waves  to  propagate  vertically  in 
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an  easterly  mean  zonal  flow;  that  is,  through  a  critical  line. 
Therefore,  when  a  propagating  wave  is  incident  upon  a  critical  line 
the  wave  energy  must  be  absorbed,  reflected,  or  possibly 
overreflected  if  the  critical  line  is  in  a  region  of  negative  qy 

(Lindzen  et  al.,  1980).  In  general,  the  vicinity  of  the  critical 
line  is  a  region  of  strong  wave  transience  and  dissipation  (Matsuno 
and  Nakamura,  1979). 

In  this  model  the  critical  region  is  generally  a  region  of 

partial  if  not  full  absorption  of  wave  energy.  This  fact  is 
demonstrated  by  the  persistence  of  the  critical  line.  The  Rayleigh 
friction  incorporated  into  the  model  tends  to  force  the  mean  flow 
back  to  its  equilibrium  state.  Since  the  critical  line  remains 

persistent  for  such  long  periods,  as  is  evidenced  in  Fig  4.3,  some 
wave  energy  must  be  absorbed  at  the  critical  level  in  order  to 

maintain  the  critical  level. 

Matsuno  and  Nakamura  (1979)  and  Schoeberl  (1980,  1981b)  have 

discussed  the  Eulerian  and  Lagrangian-mean  circulations  associated 
with  an  idealized  absorbing  critical  line  oriented  at  various  angles 
with  the  vertical.  In  Fig  4.21  are  shown  the  mass  stream  functions 
associated  with  the  two  different  formulations  as  taken  from  Matsuno 
and  Nakamura  (1979)  for  a  horizontal  critical  line.  In  the  Eulerian 
mean  formulation  there  is  rising  motion  in  the  north  with  descending 
motion  in  the  south  and  a  general  southward  transport.  This  flow 
regime  is  viewed  as  the  mean  field  response  to  the  northward  heat 
flux  of  the  planetary  waves. 
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The  Lagrangian-mean  stream  function  is  much  different  with 
singular  northward  flow  at  the  critical  line  and  weaker  return  flow 
both  above  and  below  the  critical  line.  Schoeberl  (1981b)  found 
that,  except  for  vertical  critical  lines,  there  is  always  singular 
flow  at  a  purely  absorbing  critical  line.  These  authors  use  these 
idealized  models  to  explain  stratospheric  warming  phenomena  and 
commented  on  the  consequences  of  critical  line  flows  on  the 

distribution  of  minor  constituents  in  the  stratosphere.  Schoeberl 
(1981b)  predicted  that  where  the  critical  line  changes  its 

orientation  with  the  vertical,  the  constituent  flux  will  be  highly 
divergent  and  these  regions  will  be  an  area  of  large  constituent 
density  changes. 

In  Schoeberl's  study  it  was  assumed  that  the  critical  line 
movement  was  so  slow  that  it  was  not  necessary  to  recalculate  tne 
zonal  mean  wind  field.  Nevertheless,  it  is  possible  to  calculate 
the  acceleration  of  the  mean  zonal  wind  which  was  found  to  be 
negative.  That  is,  for  the  case  of  an  absorbing  critical  line  there 
is  a  tendency  for  the  critical  line  to  descend.  Therefore  the 

transport  described  above  occurs  when  the  mean  zonal  flow  is  being 

decelerated  at  the  critical  line  and  the  critical  line  is 
descending . 

For  the  case  of  an  ascending  critical  line  the  flow  is 
accelerating  rather  than  decelerating  at  the  critical  line.  This 
being  the  case  the  sign  of  Eq .  (7)  in  Schoeberl  (1981b)  can  be 
changed  to  Indicate  that  the  critical  region  is  forcing  the  mean 
flow  to  become  more  westerly.  With  this  sign  change  a  strong 


southward  Lagrangian-mean  jet  is  found  at  the  critical  line. 

Both  the  models  of  Schoeberl  (1980,  1981b)  and  Matsuno  and 

Nakamura  (1979)  assume  the  interaction  of  the  wave  with  the  critical 
line  can  be  modeled  as  a  step  function  so  that  the  Eliassen-Palm 
(EP)  flux  is  discontinuous  across  the  critical  line*  With  such  a 
functional  dependence  of  the  wave  forcing,  the  EP  flux  divergence  is 
proportional  to  a  delta  function  and,  therefore,  the  Lagrangian-mean 
flow  is  singular  at  the  critical  line.  The  delta  function 

dependence  is  equivalent  to  putting  all  of  the  wave  transience  at 
the  critical  line.  In  the  real  atmosphere  and  in  numerical  models 

the  forcing  is  not  so  abrupt;  hence,  the  effects  of  the  critical 

line  are  spread  out  over  a  thin  layer  of  the  atmosphere. 

In  the  lagrangian-mean  framework  one  would  expect  to  see  large 
northward  fluxes  at  the  critical  line  with  the  changes  in  the  mean 
field  being  proportional  to  the  divergence  of  the  flux.  In  the 
Eulerian  formulation,  since  the  mean  velocity  is  generally 
southward,  the  strong  transport  at  the  critical  line  should  manifest 
itself  as  strong  eddy  flux  oriented  with  the  critical  line.  The 

direction  of  the  eddy  flux  will  depend  on  the  background  constituent 
density. 

In  the  following,  the  transport  for  900  m  planetary  waves  will 
be  discussed  for  both  chemical  and  conservative  transport.  It  is 
desirable  to  study  the  transport  for  conservative  tracers  because 
the  chemistry  tends  to  counter  any  changes  in  the  upper  stratosphere 
and  might  conceal  transport  mechanisms  which  would  be  important  or 
incorrect  for  other  chemical  species.  The  transport  will  be 
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discussed  with  regard  to  the  time  scale  in  Fig.  4.3  measured  from 
the  beginning  of  the  model  integration. 

The  mean  flux  is  shown  at  day  25  in  Fig  4.22.  At  all  times  the 
Eulerian  mean  flux  is  from  north  to  south  with  upward  flux  at  high 
latitudes  and  downward  flux  at  low  latitudes.  This  flux  field  is 
consistent  with  that  expected  for  a  wave  interacting  with  a  critical 
line . 

The  eddy  horizontal  flux  at  day  25  is  shown  in  Fig.  4.23  for 
conservative  transport  by  the  planetary  waves  alone.  At  this  time 
the  critical  line  has  descended  to  approximately  40  km  and  is  still 
descending.  Centered  at  the  critical  line  is  a  region  of  strong 
northward  transport  which  has  a  strong  meridional  derivative  in  both 
the  north  and  the  south.  This  region  of  northward  transport  is 
responsible  for  the  positive  anomaly  centered  at  40  km  seen  in  Fig. 
4.13.  The  region  of  southward  transport  centered  at  50  km  is 
responsible  for  the  negative  anomaly  at  50  km  in  Fig.  4.13*  The 
transport  shown  here  is  consistent  with  the  transport  as  predicted 
by  Schoeberl  (1981b)  and  Matsuno  and  Nakamura  (1979)  with  strong 
northward  transport  at  the  critical  line  and  return  flow  above. 
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Figure  4.23.  Horizontal  eddy  flux  at  day  25 
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By  day  31  (Fig  4.24a)  Che  eddy  transport  has  changed 
considerably  with  the  formation  of  several  layers  of  northward  and 
southward  flux  above  the  two  layers  that  existed  on  day  25*  The 
lowest  level  of  northward  transport  has  decreased  in  magnitude  by 
half  and  is  still  centered  at  the  critical  line  which  is  now  near 
its  maximum  descent.  At  50  km  where  there  was  strong  southward 
transport  on  day  25  there  is  now  predominately  northward 
transport.  This  reversal  of  transport  is  responsible  for  the 
decrease  in  the  magnitude  of  the  anomaly  at  50  km  around  day  30  in 
Fig.  4.13.  It  is  obvious  from  Fig.  4.13  that  the  eddy  flux  at  high 
altitudes  has  relatively  little  effect  on  the  mean  constituent 
field,  Indicating  that  the  eddy  and  mean  fluxes  are  nearly 
compensating  away  from  the  critical  line. 

In  order  to  understand  the  mechanisms  that  are  responsible  for 
the  formation  of  this  layered  transport  structure  it  is  useful  to 
consider  the  following  form  of  the  horizontal  eddy  flux  which  can  be 
derived  from  the  perturbation  continuity  equation  (Eq.  3.6)  assuming 
that  the  perturbation  experiences  exponential  growth  or  decay  and 
that  the  chemistry  is  zero. 
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This  equation  is  useful  for  diagnostics  away  from  the  critical  line. 


At  the  critical  line  many  of  the  terms  in  the  denominator  are  nearly 
zero  so  that  (4.1)  would  be  singular.  is  positive  for  a  growing 

wave . 

The  transience  at  day  31  is  shown  in  Fig.  4.24b.  Positive 
transience  is  associated  with  a  growing  wave  and  negative  transience 
is  associated  with  a  decaying  wave.  There  is  a  strong  correlation 
of  the  sign  of  the  transience  with  the  correlation  in  term  C  in  Eq . 
4.1.  When  the  wave  is  growing  (decaying)  the  correlation 
of  Re(vy)Im(w')  is  negative  (positive);  the  vertical  gradient  of 
ozone  is  negative  above  40  km,  and  the  sign  of  ku  in  the  region  of 
the  layering  is  negative.  The  horizontal  eddy  transport  in  the 
region  where  it  is  layered  is  dominated  by  term  C.  The  horizontal 
eddy  flux  in  this  region  is  proportional  to  the  vertical  constituent 
derivative.  This  dependence  is  proven  by  an  experiment  where  the 
horizontal  gradient  of  ozone  was  reversed  and  the  vertical  gradient 
held  constant.  In  this  experiment  the  layering  structure  remained 
the  same  while  the  eddy  flux  at  the  critical  line  reversed 
direction . 

The  meridional  gradient  of  the  zonal  mean  potential  vorticity, 
qy,  is  shown  in  Fig  4.24c  A  necessary  condition  for  baroclinic 
instability  is  that  this  gradient  changes  sign  (Charney  and  Stern, 
1962).  The  potential  vorticity  gradient  is  negative  in  a  large  area 
above  the  critical  line,  but  there  are  three  regions  where  it 
changes  sign  within  the  zone  of  easterlies.  The  positive 
transcience  in  Fig.  4.24b  is  localized  near  these  gradient 
changes.  Thus,  the  layered  structure  appears  to  be  due  to  the 
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transport  by  an  unstable  mode  In  the  zone  of  easterlies. 

The  positive  anomaly  In  the  ozone  field  is  persistent  until 
about  day  47  when  the  anomaly  begins  to  decrease  and  nearly 
disappears  (Figures  4.13  or  4.16).  This  feature  is  coincident  with 
the  rise  in  the  critical  line  which  occurs  around  day  50.  The 
horizontal  eddy  flux  shown  in  Fig.  4.25  indicates  that  there  is 
large  southward  flux  associated  with  the  rising  critical  line.  As 
the  critical  line  descends  again  after  day  50  strong  northward  flux 
is  once  again  associated  with  the  critical  line. 

In  actuality  the  critical  line  rises  slowly  for  two  weeks 
before  rising  rapidly  on  day  47.  During  the  time  the  critical  line 
is  nearly  stationary  and  rising  slowly,  there  is  a  small  flux  in  the 
southward  direction.  This  explains  why  there  is  not  enormous 
buildup  In  the  northern  region  while  the  critical  line  is  nearly 
steady . 

The  forcing  of  the  mean  flow  in  the  vicinity  of  the  critical 
line  is  show  In  Fig.  4.26  where  the  EP  flux  divergence  and  critical 
line  at  the  center  of  the  channel  Is  shown  for  a  dynamics  run 
similar  to  the  one  used  in  the  transport  study.  Different  dynamics 
were  used  here  merely  to  conserve  computer  resources  and  all  of  the 
essential  features  of  the  transport  dynamics  are  represented.  While 
the  critical  line  is  rising  there  is  convergence  in  the  critical 
region  which  tends  to  accelerate  the  mean  zonal  flow  in  a  westerly 
fashion.  When  the  critical  line  is  at  its  peak  on  day  42  the  EP 
flux  divergence  is  small  and  as  the  critical  line  descends  the  EP 
flux  divergence  becomes  positive  indicating  absorption  of  wave 


energy  and  a  deceleration  of  the  mean  flow. 

The  effects  of  photochemistry  are  clearly  exhibited  in  Fig. 
4.27  which  shows  the  horizontal  eddy  flux  at  day  31,  the  same  time 
as  Fig  4.24a.  The  multilayer  structure  which  was  present  in  the 
conservative  experiment  is  not  present  due  to  the  strong 

photochemistry  above  40  km.  This  indicates  that  caution  must  be 

exercised  if  the  current  dynamics  were  to  be  applied  to  other  gases 
with  different  chemistries.  If,  for  instance,  the  transition  layer 
were  located  in  the  region  where  the  unstable  mode  was  present  large 
unrepresentative  mean  field  changes  might  be  calculated. 

The  strong  poleward  transport  at  the  critical  line  is  enhanced 
by  the  photochemistry  because  the  weaker  equatorial  return  transport 
is  largely  countered  by  the  photochemistry.  The  photochemistry 

limits  the  transport  largely  to  the  conservative  region  and  the 

lower  part  of  the  transition  region  for  the  900  m  topography.  The 
Importance  of  the  critical  line  descending  below  the  transition 
region  is  demonstrated  by  the  transport  with  the  600  m  topography 
which  is  only  a  third  of  that  associated  with  the  900  m 
topography.  The  mechanism  of  transport  when  photochemistry  is 
included  is  consistent  with  the  critical  line  transport  discussed 
above.  The  high  latitude  maximum  in  Fig.  4.20  on  March  1  is 
associated  with  jet  like  transport  of  ozone  as  the  critical  line 
descends  on  day  70. 
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DAY  47.  CONSERVATIVE 


Figure  A. 25*  Eddy  horizontal  flux  at  day  A7  for  conservative 
transport  during  the  rise  of  the  critical  line  (ppm  m  sec-1), 
(bold  solid  line  is  the  critical  line) 


solid  lines  >  0 
dashed  lines  <  0 
chain-dot  “0 


Figure  A. 26.  Eliassen-Palm  (EP)  flux  divergence  during 
the  rise  of  the  critical  line,  (negative  values 
correspond  to  westerly  acceleration) 
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DAY  31.  CHEMISTRY 
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Figure  A. 27.  Horizontal  eddy  flux  on  day  31 


for  chemical  transport  (ppm  m  sec  *). 
(bold  solid  line  is  the  critical  line) 
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4. 7  Summary 

The  effect  of  the  planetary  waves  during  a  sudden  warming  is  to 
transport  ozone  poleward.  This  transport  is  explained  by  the 
mechanisms  revealed  by  Matsuno  and  Nakamura  (1979)  and  Schoeberl 
(1981b).  The  transport  by  such  strongly  transient  waves  is  much 
more  vigorous  and  confined  to  a  smaller  vertical  area  than  transport 
by  steady  state  waves. 

Northward  transport  coincides  with  a  descending  critical 
line.  Southward  transport  occurs  when  the  critical  line  ascends. 
The  magnitude  of  the  transport  varies  with  the  rate  of  ascent  or 
descent,  an  indicator  of  the  wave  absorption  properties  of  the 
critical  line.  When  the  critical  line  is  stationary  the  transport 
is  small. 

Photochemistry  enhances  the  buildup  of  ozone  near  the  pole  by 
reducing  the  return  flux  associated  with  the  Lagrangian-mean  jet  at 
the  critical  line.  The  magnitude  of  the  poleward  transport  is 
highly  dependent  on  the  critical  line  descending  from  the 
photochemical  region  into  the  transition  and  conservative  regions. 

The  diabatic  circulation  transports  ozone  northward  and 
downward  but  the  northward  transport  is  relatively  ineffective.  The 
downward  transport  produces  increases  in  ozone  in  the  lower 
stratosphere  and  replenishes  the  ozone  in  the  southern  region  that 
is  transported  northward. 

Ihe  8-plane  model  predicts  a  total  ozone  distribution  at  high 
latitudes  consistent  with  observations  but  is  less  accurate  at 


reproducing  the  vertical  structure  of  ozone.  Much  of  the  discrepancy 


can  be  attributed  to  the  lack  of  small  scale  dynamics  in  the  model 
lower  stratosphere.  In  the  next  chapter  the  effects  of  the 
simplifying  assumptions  on  the  calculated  transport  as  well  as  the 
applicability  of  the  steady  state  and  time  dependent  model  to  the 
real  atmosphere  will  be  discussed. 
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Chapter  5 


SUMMARY  AND  A  POSSIBLE  OZONE  TRANSPORT  SCENARIO 

5.1  Introduction 

In  the  preceding  chapters  o2one  transport  by  the  diabatic 
circulation  and  both  steady  state  and  time  dependent  planetary 
waves  on  a  &-plane  has  been  described.  The  results  simulate  many  of 
the  large  scale  features  of  the  wintertime  ozone  distribution  with 
surprising  accuracy.  In  this  chapter  the  effects  of  some  of  the 

simplifying  assumptions  incorporated  into  the  model  will  be 
evaluated,  and  the  results  will  be  Interpreted  with  regard  to  their 
applicability  to  the  atmosphere.  It  should  be  remembered  that  the 
purpose  of  the  study  Is  primarily  to  Investigate  transport 

mechanisms,  not  to  accurately  model  all  aspects  of  the  ozone 
distribution . 

Many  simplifications  were  made  both  tacitly  and  explicitly  in 
order  to  facilitate  both  the  execution  and  the  interpretation  of  the 
models.  These  specifically  include  the  omission  of  small  scale 
dynamics  in  the  lower  stratosphere;  thi  lack  of  temporal  and 

realistic  latitudinal  variation  in  the  photochemistry;  the 
specification,  rather  than  the  calculation,  of  the  diabatic 
circulation;  and  the  many  assumptions  incorporated  in  applying 
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quasigeostrophic  planetrary  waves  on  a  g-plane  to  a  spherical 
atmosphere.  The  effects  of  these  assumptions  on  the  final  results 
are  varied  and  difficult  to  evaluate  quantitatively  but  in  some 
instances  their  qualitative  influences  can  be  described. 

5.2  Low  latitude  ozone 

One  of  the  most  serious  inadequacies  of  the  model  is  the 
overestimation  of  planetary  wave  transport  out  of  the  southern 
region  of  the  model  during  sudden  warmings.  This  excessive 
transport  is  exacerbated  by  the  underestimation  of  the  total  ozone 
south  of  50°  on  Dec.  15,  the  time  when  the  planetary  waves  are 
turned  on. 

a.  The  diabatic  circulation  and  the  initial  ozone  distribution 

The  problem  of  too  little  ozone  on  Dec.  15  might  be  attributed 
to  the  Initial  ozone  distribution  being  in  error  either  with  regard 
to  total  amount  or  the  specification  of  mean  gradients  (see  Section 
4. 2d).  If  the  initial  distribution  is  assumed  to  be  correct,  then 
the  underestimation  of  ozone  on  Dec.  15  might  be  due  to  the  improper 
designation  of  either  the  horizontal  or  the  vertical  diabatic 
velocity  in  the  southern  region.  If  the  vertical  velocity  is  too 
small,  the  transport  from  high  altitude  source  regions  to  low 
altitudes  would  be  underpredicted. 

At  the  southern  boundary  there  is  a  northward  flux  of  ozone 
associated  with  the  meridional  component  of  the  diabatic  velocity 
during  fall  and  winter.  If  this  component  is  too  small,  then  the 
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flux  across  Che  southern  boundary  may  be  underestimated. 
Furthermore,  the  ozone  has  been  assumed  to  be  constant  at  the 
southern  boundary  point  (37.5°);  when  in  fact,  there  is  a  102 
increase  of  total  ozone  mixing  ratio  during  the  winter  (see  Fig. 
2-8).  The  specification  of  a  constant  mixing  ratio  at  the  southern 
boundary  is,  therefore,  not  strictly  correct. 

The  inclusion  of  planetary  waves  in  forming  the  Dec.  15  ozone 
field  would  tend  to  increase  the  total  ozone  error  at  low  latitudes 
since  both  steady  state  and  transient  waves  cause  a  depletion  of 
ozone  in  the  southern  part  of  the  8-channel .  Even  with  an  improved 
estimate  of  the  ozone  field,  the  transport  of  ozone  out  of  the  south 
during  the  initial  warming  is  too  strong.  There  is  a  large  poleward 
transport  out  of  the  southern  region  associated  with  the  initial 
warming  which  causes  a  more  than  20  DU  decrease  in  total  ozone.  Any 
decrease  of  ozone  at  40°  N  during  winter  is  unrealistic. 

b.  Geometric  effects 

The  geometric  effects  of  the  8-plane  would  also  contribute  to 
the  overestimation  of  transport  out  of  the  south.  Figure  5.1  shows 
a  spherical  atmosphere  divided  into  cells  and  the  corresponding 
cells  on  a  8-plane.  The  high  latitude  cells  on  the  sphere  contain 
less  volume  than  the  low  latitude  cells;  in  fact,  a  latitude  band 
2.15°  wide  (the  approximate  horizontal  grid  length)  at  82.5° 
contains  less  than  a  third  the  volume  of  a  similar  cell  centered  at 
60°.  Likewise,  a  cell  at  60°  contains  only  64%  of  the  volume  of  a 
cell  at  37.5°.  Since  the  8-plane  geometry  assumes  a  constant  cell 
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volume,  the  constituent  divergence  at  low  latitudes  is  overestimated 
while  the  constituent  convergence  at  high  latitudes  is 
underestimated.  In  other  words,  the  ozone  transported  north  should 
be  distributed  over  a  smaller  volume  than  it  is,  and  the  ozone 
removed  from  the  south  should  be  drawn  from  a  larger  volume  than  it 
is. 

In  order  to  understand  the  nature  of  this  error  consider  the 
meridional  term  of  the  mean  constituent  divergence 


rise  h  (^C088) 
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(5.1) 


In  the  3-plane  approximation  cos0  is  replaced  by  cos0o,  where  0O  is 
the  latitude  of  the  center  of  the  channel,  and  (5*1)  becomes 

rise  le  <^c°s0>  -  (y-a0)  (5.2) 

So  it  is  seen  that  the  3-plane  estimate  is  exact  only  at  the  center 
of  the  channel.  The  error  due  to  the  geometric  effect  is  related  to 
the  second  term  on  the  right  hand  side  of  (5.1),  and  is  proportional 
to  the  tangent  of  the  latitude  and  the  magnitude  of  the  flux.  At 
very  high  latitudes  the  tangent  term  would  become  very  large.  In 
this  model  both  the  mean  and  eddy  horizontal  fluxes  approach  zero  at 
the  northern  boundary  (82.5°);  so  that,  the  errors  remain  small.  At 
low  latitudes,  the  tangent  term  is  smaller  than  that  at  the  channel 
center;  therefore,  the  divergence  would  be  overestimated  on 
the  0-plane.  The  fluxes  associated  with  the  planetary  wave  dynamics 
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are  small  at  the  southern  boundary  because  the  wave  amplitude  is 
assumed  to  be  zero.  However,  the  flux  associated  with  the  dlabatlc 
circulation  at  the  southern  boundary  Is  relatively  large  and  could 
contribute  to  the  underestimation  of  the  ozone  in  the  south. 

c.  Photochemistry 

The  photochemistry  is  not  truly  conservative  in  the  lower 
stratosphere,  and  sufficiently  rapid  photochemistry  might  exist  at 
southern  latitudes  to  affect  the  transport.  Figure  5.2  from  Gille 
et  al .  (1980)  shows  an  observational  estimate  of  the  transition 

region  based  on  the  phase  between  the  geopotential  and  the  ozone 
wave  (Hartmann  and  Garcia,  1979).  In  the  southern  parts  of  this 
figure  the  transition  region  descends  to  nearly  20  km.  Inclusion  of 
chemistry  between  28  and  30  km  would  help  to  alleviate  the  severe 
reduction  calculated  in  the  southern  part  of  the  model  by  generating 
new  ozone  to  replace  that  transported  northward . 
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Figure  5.2.  Latitudinal  dependence  of  the  lower  boundary 
of  the  photochemical ly  controlled  region  (circles) 
and  the  upper  boundary  of  the  dynamically 
controlled  region  (squares),  for  November  and 
December  (from  Gille  et  al.,  1980). 
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Figure  5.3.  Meridional  cross  section  of  the  zonal 

wind  velocity  (m  sec-*),  (from  Palmer,  1981) 


d.  Planetary  waves 

The  excessive  transport  out  of  the  southern  regions  of  the 
model  during  sudden  warmings  is  influenced  by  all  of  the  processes 
discussed  above;  however,  the  problem  is  exaggerated  by  the 
simplicity  of  the  planetary  wave  dynamics.  The  critical  lines  in 
many  of  the  figures  in  Chapter  4  (Fig  4.23  for  instance)  are  nearly 
horizontal  and  extend  from  the  northern  boundary  to  the  southern 
boundary.  During  stratospheric  sudden  warmings,  the  critical  line 
generally  does  not  extend  south  of  50°,  which  is  close  to  the  center 
of  this  model  (Fig.  5.3).  Furthermore,  the  critical  line  is  rarely 
perfectly  horizontal  and  is  likely  to  be  in  the  photochemical  region 
at  midlatitudes  and  in  the  conservative  region  at  high  latitudes. 
Therefore,  much  of  the  excessive  transport  out  of  the  south  exists 
because  the  model  distorts  the  horizontal  scale  of  the  warming 
event.  More  realistic  initial  wind  profiles  would  help  alleviate 
this  problem. 

5.3.  Vertical  structure 
a.  Small  scale  dynamics 

In  Chapter  4  the  high  latitude  total  ozone  was  well  estimated, 
but  the  polar  buildup  occurred  higher  than  observed.  This  was 
attributed  to  the  lack  of  subplanetary  scale  dynamics  in  the  lower 
stratosphere.  Since  ozone  is  nearly  conservative  in  the  lower 
stratosphere,  transport  here  would  serve  only  to  redistribute  ozone 
and  not  change  the  total  amount.  Planetary  waves  which  are  not 
capable  of  propagating  to  very  high  altitudes  (wavenumber  2-5)  are 
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often  present  in  the  lower  stratosphere  as  are  the  remnants  of 
tropospheric  weather  systems.  Also,  it  is  well  known  that  intense 
baroclinic  zones  in  the  upper  troposphere  cause  excursions  of  ozone 
into  the  troposphere  and  may  be  instrumental  in  the  downward 
transport  of  ozone  (Kida,  1977;  Shapiro,  1980).  Thus,  it  appears 
that  a  variety  of  mechanisms  °xist  which  can  redistribute  ozone 
downward,  while  maintaining  the  column  density,  and  that  the 
inability  of  the  current  model  to  properly  simulate  the  vertical 
distribution  does  not  vitiate  the  importance  of  the  observed 
transport  mechanisms. 

It  might  be  possible  to  model  the  subplanetary  processes  in  the 
lower  stratosphere  related  to  the  redistribution  of  ozone  by 
diffusion.  However,  for  chemically  active  constituents  other  than 
ozone,  many  of  the  problems  that  appear  when  applying  diffusion 
models  to  ozone  might  become  important.  In  order  to  evaluate  that 
effects  of  pollutants  on  ozone,  the  dynamics  of  the  pollutants  must 
be  accurately  calculated;  therefore,  any  diffusion  model  must  be 
used  judiciously. 

b.  Photochemistry 

Many  aspects  of  the  chemistry  of  ozone  as  used  In  this  model 
are  inadequate  and  might  contribute  to  incorrect  results.  The 
chemistry  was  deliberately  simplified  so  that  attention  could  be 
focused  on  the  transport  and  interaction  of  the  dynamics  and 
photochemistry.  Rather  than  speculate  on  the  particular  effects  of 
the  omission  of  certain  chemical  species,  it  is  more  pertinent  to 
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comment  on  the  spatial  variation  of  the  simple  chemistry  included. 

Garcia  and  Hartmann  (1980)  found  that  the  latitudinal  variation 
of  photochemistry  (the  effect  of  the  polar  night)  was  important  in 
enhancing  the  polar  buildup  of  ozone.  In  Section  4.3  the 
latitudinal  variation  of  photochemistry  was  shown  to  have  little 
effect  on  transport  by  the  diabatic  circulation  into  the  lower 
stratosphere.  When  ozone  is  assumed  to  be  conserved  in  the  polar 
night,  a  significant  depletion  develops  over  the  course  of  the 
winter  above  40  km.  Very  weak  chemistry,  though,  will  significantly 
moderate  this  depletion  since  the  advection  time  scales  are  so  long. 

As  the  planetary  wave  transport  is  modeled  here  for  time 
dependent  waves,  the  latitudinal  variation  of  photochemistry  would 
have  relatively  little  effect  at  high  latitudes.  Most  of  the 
transport  occurs  below  the  transition  layer  where  ozone  is  nearly 
conservative  at  all  times.  The  reduction  of  chemistry  would  have  a 
more  noticeable  effect  on  transport  during  minor  warmings  when  the 
critical  line  does  not  descend  very  low. 

5.4  Unstable  waves 

During  the  transient  transport  an  unstable  wave  exists  within 
the  critical  level  which  causes  the  layered  transport  phenomemon 
discussed  in  Section  4.6.  For  ozone  this  layering  had  little 
effect,  because  it  occurred  in  the  photochemical  region.  For  a 
conservative  tracer  or  a  compound  with  its  transition  region  located 
in  the  unstable  region,  the  unstable  wave  may  have  considerable 
transport  consequences. 
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5.5  Summary  and  discussion 

Despite  the  simplifications  discussed  above,  the  transport  as 
calculated  in  this  model  reproduces  many  of  the  features  of  the 
wintertime  ozone  distribution.  Furthermore,  due  to  the  simplicity 
of  the  model,  specific  mechanisms  which  are  responsible  for 
transport  can  be  identified. 

a.  Steady  state  waves 

In  Chapter  3,  transport  by  steady  state  dissipative  planetary 
waves  was  examined  from  both  an  Eulerian  and  Lagrangian-mean 
viewpoint.  In  the  lower  stratosphere  these  waves  were  found  to 
cause  ozone  depletion  at  low  latitudes  and  an  ozone  enhancement 
toward  the  pole.  The  increase  in  ozone  in  the  northern  region  is 
consistent  with  the  observations;  however,  there  is  not  an  observed 
decrease  at  low  latitudes  as  is  calculated  here.  The  ozone 
transport  in  the  lower  stratosphere  was  due  to  advection  by  the 
Lagrangian-mean  velocity  associated  with  the  planetary  waves. 

In  the  lower  part  of  the  transition  region  there  is  a  increase 
of  mean  ozone  at  high  latitudes  and  a  decrease  at  low  latitudes  (see 
Figures  3.7  and  3.14).  This  behavior  in  the  transition  region  was 
shown  to  be  due  to  the  strong  photochemical  gradients  in  the  region 
(Sec.  3.3d).  The  particle  trajectories  associated  with  the 
planetary  waves  are  such  that  parcels  move  through  a  region  of 
widely  varying  photochemistry.  As  a  result  the  parcels  experience  a 
photochemical  forcing  different  from  the  mean  field  chemistry  at  the 
parcel's  equilibrium  position.  Thus  the  perturbation  chemistry 


cannot  be  balanced  by  the  mean  field  chemistry,  and  changes  In  the 
mean  ozone  mixing  ratio  occur. 

Where  the  photochemical  gradients  are  small  there  is  little  or 
no  effect  of  the  planetary  waves  on  the  mean  field.  The  steady 
state  model  indicates  that  transport  by  steady  planetary  waves  might 
be  adequately  represented  by  advectlon  by  the  Lagrangian-mean 
velocity  in  the  lower  stratosphere,  and  some  scheme  which  Includes 
the  parcels  orbiting  through  regions  of  strong  photochemical 
gradients.  Indeed,  a  model  based  on  observed  trajectories  has  been 
recently  proposed  by  Garcia  and  Solomom  (1982). 

An  accurate  estimate  of  the  Lagrangian-mean  velocity  for  a 
realistic  model  of  the  stratosphere  is  difficult  to  calculate  since 
the  Lagrangian-mean  velocity  Is  the  residual  between  two  large 
terms.  In  atmospheric  data  since  there  are  large  uncertainties  in 
the  estimates  of  the  quantities  needed  to  produce  the  Lagrangian- 
mean  displacement  fields,  it  would  be  even  more  difficult  to  form 
meaningful  estimates  of  the  Lagrangian-mean  velocity  from 
observations . 

In  Section  3.3c  the  steady  state  ozone  transport  using  the 
Eulerian  residual  circulation  was  calculated.  The  results  were 
similar  to  those  found  using  the  Lagrangian-mean  velocity,  but  the 
advectlon  was  consistently  overestimated.  The  greatest  differences 
between  the  residual  velocity  and  the  Lagrangian-mean  velocity  occur 
in  the  photochemical  region  so  that  the  chemistry  obscures  the 
errors  in  transport.  For  a  constituent,  with  a  different 
photochemistry  from  ozone,  there  might  be  significant  differences 
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between  the  Lagrangian-mean  and  residual  transport.  It  appears  that 
the  Eulerian  formulation  is  best  for  carrying  out  transport 
calculations  in  realistic  models,  and  the  other  formulations  are 
useful  for  Interpreting  the  results. 

Over  the  course  of  the  winter,  the  computed  ozone  increases  in 
the  lower  stratosphere  at  high  latitudes  show  reasonable  agreement 
with  observation.  Using  the  "realistic"  zonal  wind  profile  of 
Schoeberl  (1981a)  too  much  ozone  was  transported  into  the  northern 
region.  The  depletion  of  ozone  in  the  southern  lower  stratosphere, 
could  be  partially  alleviated  with  the  inclusion  of  the  diabatic 
circulation.  The  diabatic  circulation  would  serve  to  transport 
ozone  into  the  lower  stratosphere  and  help  replenish  the  depleted 
areas.  However,  even  with  the  inclusion  of  the  diabatic  circulation 
the  steady  state  ozone  model  would  not  describe  the  winter  ozone 
transport  and  would  exhibit  many  of  the  failures  of  the  total  ozone 
field  calculated  with  the  diabatic  circulation  alone  (Sec.  4.3). 
The  north-south  transport  associated  with  the  planetary  waves  is 
larger  than  that  associated  with  the  diabatic  circulation  but  is 
still  not  adequate  to  supply  the  required  northward  transport 
observed  in  the  winter. 
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Figure  5.4.  Simulated  zonal  mean  total  ozone  (DU) 
(from  Mahlman  et  al . ,  1980). 
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It  is  interesting  to  consider  the  results  of  Mahlman  et  al « 
(1980)  within  the  context  of  the  simple  transport  models  used 
here.  Figure  5.4  taken  from  Mahlman  et  al .  shows  their  predicted 
total  ozone  field.  Comparison  with  Fig.  2.8  shows  that  they  predict 
a  spring  ozone  maximum  which  occurs  too  far  south,  too  late  in  the 
year,  and  that  is  too  large.  There  is  an  intrusion  of  ozone  into 
the  polar  region  during  the  springtime  reversal  of  the  mean 
circulation.  These  errors  are  consistent  with  the  errors 
associated  with  transport  by  the  diabatic  circulation  alone  or  the 
diabatic  circulation  with  steady  state  planetary  waves.  The  results 
from  this  study  suggest  that  it  is  the  inability  of  the  Mahlman  et 
al .  model  to  simulate  stratospheric  warmings  which  is  responsible 
for  much  of  the  error  in  their  total  ozone  field. 

b.  Time  dependent  waves 

The  transport  by  transient  planetary  waves  (Chapter  4)  was 
found  to  be  much  different  from  steady  state  transport,  but  it  is 
conceptually  quite  simple.  Intense  poleward  transport  is  associated 
with  the  deceleration  and  reversal  of  the  zonal  mean  flow  during 
warmings,  and  equatorward  transport  is  associated  with  the 
reacceleration  of  the  flow  after  the  warming.  The  transport  by 
highly  transient  planetary  waves  is  largely  confined  to  the 
horizontal  plane  with  only  a  small  vertical  component.  More 
realistic  zonal  wind  profiles  might  generate  critical  lines  which 
have  a  larger  vertical  transport  associated  with  them  (Schoeberl 
1981a,  1981b). 


141 


With  the  transport  being  so  much  difierent  for  transient  and 

steady  state  waves,  the  question  arises  as  to  the  consistency 
between  the  two  models.  To  compare  the  two  models  an  experiment  was 

run  in  which  the  time  dependent  model  was  forced  with  a  200  n 

mountain.  Such  forcing  causes  large  wave-mean  flow  interaction 
initially  but  there  is  never  a  reversal  in  the  zonal  mean  wind 

field,  and  an  equilibrium  state  is  reached  with  westerlies  present 
in  the  entire  domain.  At  an  altitude  of  40  km  the  wave  amplitude 
assumes  a  near  steady  value  of  approximately  350  gpm. 

The  anomaly  field  is  shown  at  74°  in  Fig.  5.5.  There  is  a 
positive  change  in  the  lower  part  of  the  transition  region  and  a 

negative  change  just  below  the  transition  region  which  is  in 

qualitative  aggreement  with  the  steady  state  model  (Fig  3.14).  The 
buildup  found  in  the  lower  stratosphere  at  high  latitudes  is  not 
reproduced  in  this  experiment.  The  changes  in  the  mean  field  for 
the  200  m  waves  are  an  order  of  magnitude  or  more  smaller  than  those 
calculated  with  the  600  and  900  m  waves  and  are  also  much  smaller 
than  those  calculated  with  the  steady  state  model. 

The  difference  in  the  lower  stratosphere  between  the  steady 
state  and  time  dependent  models  is  due  to  wave-mean  flow  interaction 
in  the  upper  stratosphere.  In  the  time  dependent  model  the  wave 
induced  mean  velocity  fields  are  strongest  in  the  region  of  large 
wave-mean  flow  interaction.  These  fields  are  active  at  a 

significantly  higher  altitude  than  in  the  steady  state  model  so  that 

the  largest  advection  by  the  planetary  waves  does  not  occur  in  the 
lower  stratosphere.  This  implies  that  the  assumption  in  the  steady 


state  model  (Section  2.3a)  that  the  changes  in  the  mean  field  are 
small  may  significantly  affect  the  results. 

It  is  interesting  to  compare  the  eddy  flux  field  from  the  200  m 
experiment  with  those  actually  observed  by  Gille  et  al .  (1980)  (Fig 
5.6).  The  model  results  agree  well  with  the  observations  suggesting 
that  the  200  m  experiment  might  be  appropriate  for  the  early  winter. 

The  steady  state  model  may  be  useful  for  modeling  transport  in 
the  late  fall  and  early  winter;  however,  the  time  dependent  model 
more  accurately  represents  the  atmospheric  transport  in  late  winter 
and  spring.  Figure  5.7  schematically  summarizes  the  different 
transport  mechanisms  studied  here.  Both  planetary  wave  models 
deplete  ozone  in  the  south  and  enhance  it  in  the  north.  The  effects 
of  the  transition  layer  are  not  obvious  for  transport  during 
stratospheric  warmings  as  is  expected  because  the  transport  at  the 
critical  line  Is  so  strong  that  it  obscures  this  relatively  weak 
effect . 

The  different  topographies  used  to  generate  the  time  dependent 
dynamics  in  Chapter  4  allow  for  the  study  of  transport  during  all 
types  of  warming  events:  major  (900  m) ,  minor  (600  m),  and  final. 
The  critical  line  remains  above  32  km  in  the  600  m  experiment  and  is 
analogous  to  a  minor  warming  while  the  900  m  critical  line  extends 
below  32  km  as  it  does  during  a  major  warming.  The  persistence  of 
the  critical  line  for  such  a  long  period  of  time  after  the  900  m 
warming  allows  for  the  study  of  transport  during  a  final  warming. 
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b.  diabattc  circulation  and  transient  waves 

Figure  5.7.  Schematic  summary  of  transport  mechanisms. 
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The  transport  for  the  900  m  topography  is  a  factor  of  3  or  more 
greater  than  that  for  the  600  m  topography.  The  increased  transport 
associated  with  the  900  m  topography  arises  because  the  critical 
line  descends  through  the  transition  region  into  the  conservative 
region  where  the  transport  is  not  countered  by  the  chemistry.  This 
suggests  that  transport  during  major  warmings  should  be  more 
pronounced  than  during  minor  warmings.  This  behavior  is  consistent 
with  the  observations  in  Fig.  2.10  where  in  the  two  years  that  major 
warmings  occur  the  ozone  maximum  at  high  latitudes  is  6%  greater 
than  the  year  in  which  no  major  warming  occurs. 

After  both  major  and  minor  warmings  there  is  a  tendency  for 
westerlies  to  develop  and  the  stratosphere  to  return  to  its 
climatological  state.  In  this  event  the  critical  line  rises  and 
much  of  the  ozone  transported  into  the  north  is  transferred  back  to 
the  south.  This  is  clearly  seen  by  the  slight  rise  in  the  critical 
line  that  occurs  in  the  900  ra  experiment  on  day  47  (Fig.  4.3).  Such 
behavior  might  be  responsible  for  the  oscillatory  structure  seen  in 
Fig.  2.10. 

The  rise  in  the  900  m  critical  line  (Fig  4.3)  followed  by  the 
descent  and  persistence  of  the  critical  line  is  representative  of 
what  happens  in  a  winter  when  there  is  a  major  warming  followed  by  a 
final  warming.  The  descent  of  the  critical  line  after  Feb.  1  is 
Instrumental  in  the  formation  of  the  sharp  high  latitude  total  ozone 
maximum  calculated  with  the  combined  circulation  (Fig  4.9). 

The  suggestion  is  made  from  this  model  that  transport  during 
major  sudden  warmings  is  large  and  may  have  significant  effect  on 
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the  zonal  mean  ozone  distribution,  but  that  the  high  latitude  spring 
maximum  is  due  largely  to  transport  during  the  spring  reversal. 
Following  the  spring  reversal  the  stratosphere  remains  in  an 
easterly  state  so  that  there  is  not  a  strong  reverse  flow.  It  might 
also  be  suggested  that  the  fall  minimum  is  associated  with  the 
reformation  of  westerlies  in  the  stratosphere. 

c.  A  possible  transport  scenario 

Using  the  transport  mechanisms  described  in  this  dissertation, 
and  the  observed  behavior  of  ozone  in  both  the  stratosphere  and 
general  circulation  models;  it  is  possible  to  derive  a  possible 
scenario  for  large  scale  ozone  transport  in  the  stratosphere.  First 
it  is  necessary  to  form  an  analogue  of  the  dynamics  in  the  different 
models  used  here  to  the  atmosphere. 

Diabatic  forcing  determines  the  direction  of  the  mean  zonal 
wind,  which  in  turn  controls  the  propagation  of  waves  into  the 
stratosphere;  therefore,  the  planetary  wave  transport  and  the 
diabatic  transport  are  not  completely  independent  as  modeled  here. 
The  changeover  from  winter  westerlies  to  summer  easterlies  after  the 
final  stratospheric  warming  is  related  the  insolation  at  high 
latitudes  forcing  the  zonal  mean  wind  in  a  easterly  direction.  It 
is  possible  that  major  and  minor  warmings  in  the  winter  cause  the 
springtime  reversal  to  occur  somewhat  earlier  than  if  diabatic 
forcing  alone  were  responsible  for  the  reversal.  This  possibility 
is  supported  by  the  reversal  of  the  zonal  westerlies  in  the  Southern 
Hemisphere  which  occur  somewhat  later  in  the  winter  than  in  the 
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Northern  Hemisphere.  Stratospheric  warmings  are  not  as  a 
predominate  feature  of  the  Southern  Hemisphere  circulation  as  of  the 
Northern  Hemisphere  circulation.  Furthermore,  the  model  of  Mahlman 
et  al .  (1980)  shows  the  intrusion  of  ozone  into  high  latitudes 

coincident  with  the  spring  reversal.  This  intrusion  occurs  later 
than  is  actually  observed.  Since  this  general  circulation  model 
does  not  simulate  stratospheric  warmings,  the  idea  that  planetary 
wave  activity  is  important  in  the  formation  of  the  summertime 
easterlies  is  further  supported. 

During  summer,  with  continual  dlabatic  forcing  at  the  pole, 
easterlies  dominant  the  stratospheric  circulation  and  planetary 
waves  become  unimportant.  In  the  fall  as  the  insolation  decreases 
at  high  latitudes  westerlies  reform  and  planetary  wave  activity  once 
again  increases. 

In  Figure  4.8  the  diabatic  circulation  produced  a  fall  minimum 
in  October  in  good  agreement  with  observation.  The  minimum  occurs 
at  the  time  that  diabatic  circulation  is  changing  over  from  upward 
motion  at  high  latitudes  during  the  summer  to  the  downward  motion 
present  in  the  winter.  This  minimum  is  also  coincident  with  the 
change  over  from  zonal  mean  easterlies  to  zonal  mean  westerlies  in 
the  stratosphere. 

After  the  ozone  minimum  in  October,  downward  advection  by  the 
diabatic  circulation  tends  to  build  up  ozone  in  the  lower 
stratosphere.  As  ozone  is  increasing,  planetary  wave  activity  is 
becoming  more  vigorous  and  warming  events  begin  to  occur.  "Steady" 
waves  constantly  transport  ozone  poleward,  but  during  major 
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stratospheric  warnings  large  amounts  of  ozone  are  transported 
northward.  Due  to  small  scale  "diffusive"  processes,  some  of  the 
ozone  transported  to  high  latitudes  during  warmings  remains  there 
leading  to  the  enhanced  extratropical  total  ozone  evidenced  in  Fig. 
2.10.  Also,  as  ozone  is  being  advected  into  the  conservative 
region,  small  scale  dynamics  conservatively  redistribute  it  into  the 
lower  stratosphere  and  troposphere.  All  of  these  features  are 
evident  in  Fig.  2.6  where  there  is  first  an  increase  of  ozone  at  30 
mb  and  then  downward  propagation  from  mid-January  through  March. 

At  the  time  of  the  final  warming,  which  is  determined  not  only 
by  planetary  wave  forcing  but  also  by  Insolation  at  high  latitudes, 
there  is  strong  northward  transport.  The  total  o^one  maximum  is 
coincident  with  the  breakdown  of  the  polar  vortex;  this  maximum 
occurs  somewhat  earlier  than  the  maximum  calculated  using  t-’ne 
diabatic  circulation  in  Fig.  A. 9.  The  transport  during  the  final 
warming  is  strong  because  of  the  interaction  of  the  well  developed 
winter  planetary  waves  with  the  critical  line.  During  fall  when  the 
westerlies  reappear,  such  dramatic  transport  does  not  occur  because 
the  rave-mean  flow  interaction  is  much  smaller.  After  the  spring 
maximum,  during  summer,  the  diabatic  circulation  slowly  reduces  the 
total  ozone  as  it  transports  it  upward  into  the  photochemical  region 
where  the  ozone  then  assumes  its  equilibrium  value. 

In  the  Southern  Hemisphere  the  intrusion  of  ozone  into  polar 
latitudes  occurs  during  the  reversal  of  the  mean  circulation 
associated  with  the  diabatic  forcing.  This  indicates  the  reduced 


importance  of  major  stratospheric  warmings  in  Southern  Hemisphere 
dynamics . 

Transport  as  described  by  this  scenario  suggests  that  the 
diabatic  circulation  ( Brewer-Dobson  circulation)  is  responsible  for 
many  of  the  large  scale  features  of  the  observed  ozone 
distribution.  However,  the  diabatic  circulation  cannot  effectively 
transport  ozone  into  polar  regions,  and  transport  to  very  high 
latitudes  is  accomplished  by  planetary  waves  during  sudden  and  final 
warmings.  Subplanetary  scale  and  evanescent  planetary  scales 
redistribute  ozone  in  the  lower  stratosphere  and  transport  ozone 
into  the  troposphere. 

Finally,  some  comment  must  be  made  about  the  implications  of 
the  modeled  transport  on  diffusion  models.  The  transport  by  the 
diabatic  circulation  is  advective  and  simple  to  model.  Hie 
advective  transport  by  steady  planetary  waves  does  not  seem  to  be 
adequate  to  model  the  northward  transport  of  ozone  during  winter. 
The  transport  during  sudden  warmings  is  large,  coherent,  rapid,  and 
not  diffusive. 

The  model  results  indicate  that  the  large  scale  ozone  transport 
is  dominated  by  vertical  diabatic  transport  and  rapid  pulses  of 
horizontal  transport  during  warmings.  The  planetary  waves  are  not 
responsible  for  much  vertical  transport;  though,  subplanetary  waves 
in  the  lower  stratosphere  might  cause  significant  vertical 
transport . 

The  diffusion  coefficients  as  calculated  by  Strobel  (1981) 
clearly  indicate  the  complexities  involved  in  accurately  calculating 
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diffusion  coefficients.  They  are  a  strong  function  of  source  terms, 
mean  field  gradients,  and  wave  transience.  The  type  of  transport 
calculated  here  is  highly  transient  and  indicates  that  rapidly 
varying  diffusion  coefficients  are  required  to  properly  simulate 
transport  processes.  If  these  coefficients  are  to  be  calculated 
from  dynamic  and  chemical  quantities  then  the  entire  purpose  of 
diffusion  parameterizations ,  computational  simplicity,  is  defeated. 
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Appendix  A 


NOMENCLATURE  AND  CONVENTIONS 


(u,v,w) 


Newtonian  cooling  and  Rayleigh  friction  coefficients 

Coriolis  parameter  (lx  10~^sec*‘*) 

scale  height  (7  km) 

zonal  wave  number  (6.28  x  10 cm-1) 

Brunt-Vaisala  frequency  squared  (4  x  10'^sec-^) 
pressure;  pQ  reference  pressure 

horizontal  derivative  of  the  mean  potential  vorticity 

photochemical  source  term 

-  (uj,  U2»  U3)  ■  (zonal,  meridional,  vertical) 

velocity 
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(v*,  v*)  *  (meridional,  vertical)  residual  velocity 
(Andrews  and  McIntyre,  1976) 


t  time 

(x,y,z)  »  (xp  x2,  X3)  -  (east-west,  north-south, 

z  -  -Hln  —  )  coordinates 
po 

d  meridional  variation  of  Coriolis  parameter,  f 

A  photochemical  time  constant 

p  ozone  mixing  ratio 

(S.n.O  -  (S1.S2.53)  *  (longitudinal,  meridional, 

vertical) 

displacement  fields 


4  geopotential 

CONVENTIONS : 

For  a  given  T 

?  i  zonal  average  (Euler ian) 

T'  -  T  -  f  5  zonal  deviation 

T**  =  Lagrangian-mean  (Andrews  and  McIntyre,  19  78) 
T8  •  T**  -  T  £  Stokes  correction  (see  Eq .  (3.5)) 
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Appendix  B 


NUMERICAL  PROPERTIES  OF  THE  ADVECTION-DIFFUSION  EQUATION 


The  one  dimensional  advection-dif fusion  equation  is 


3S  as 
at  +  =  D 


a2s 

,  2 
3x 


(B.l) 


Using  second  order,  centered  time  differences  and  generalized 
spatial  differences  (B.l)  can  be  written  in  discrete  form  as 


_n+l  n-1 
bt  "  Si  ~ 


+ 


2D  At 


Ax 


-TT* < Ai(S"+r  sh> +  s;.2» 

{  A2 (Si+r  2S1+  Si-1>  +  82(SJ+2-  2Sl+  Sl-2» 


(B.2) 


where  Ax  is  the  spatial  grid  increment  and  At  is  the  time  step. 
The  superscipts  are  time  indices  and  the  subscripts  are  space 
indices.  Both  second  order  and  fourth  order  differences  will  be 
considered,  and  the  values  of  the  A  and  B  coefficients  in  (B.2)  are 
given  in  table  B.l. 


Table  B.l 


Coefficients  for  Spatial  Difference  Scheme 

2n<*  order  4t*1  order 


Ax  1/2  2/3 

0  -1/12 

A2  1  4/3 

B2  0  -1/12 


Assuming  the  solution  can  be  written  as 


S  -  BnAteimkAx  (B. 3) 

where  m  is  the  spatial  index,  then  (B.2)  can  be  written  following 
Haltiner  and  Williams  (1980,  Chapter  5)  as 
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(B.4) 


BAt  -  B_At-  1  — -AAt-  {Ajsin(kAx)  +  BlCos(2kAx)} 

+  At  {A2(cos(kAx)-l)  +  B2(cos(2kAx)-l)} 

Ax^ 

In  the  transport  model  there  is  no  diffusion  so  first  consider  (B.4) 
with  D  «  0-  It  can  then  be  shown  that  for  second  order  differences 

<  1  for  stability,  and  for  fourth  order  differences 

<  *73  for  stabilty. 

Thus,  for  a  given  Ax,  the  fourth  order  scheme  requires  a  smaller 
time  step  than  the  second  order  scheme  to  assure  numerical 
stability.  However,  as  demonstrated  in  Fig.  2.1,  the  phase  error 
for  the  fourth  order  differences  is  considerably  reduced  (see 
Haltiner  and  Williams,  1980,  Table  3.1).  Haltiner  and  Williams 
(1980)  discuss  in  detail  the  numerical  properties  of  the  second 
order  and  fourth  order  advection  equation.  Krelss  and  Oliger  (1972) 
argue  that,  by  comparing  the  error  in  the  phase  speed  with  the 
number  of  computations,  the  fourth  order  scheme  is  "optimal." 

Clancy  (1981)  suggests  using  forward  time  differences  and 
second  order  spatial  differences  to  Integrate  the  advection- 
dif fusion  equation.  Clancy  derives  a  new  stability  criterion. 
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TR  <  1 
c 


where 

r  u  Ax 
At 


and 


R 

c 


u  Ax 
~D 


(B.5) 


Using  this  stability  criterion  Rc  can  be  infinitely  large  if  T  is 
small  enough  to  keep  the  inequality  (B.5)  true.  It  had  been  shown 
previously  that  Rc  must  be  less  than  or  equal  to  2,  if  this  scheme 
is  to  be  numerically  stable.  In  the  derivation  of  the  stability 
criterion  (B.5),  however,  Clancy  effectively  divided  by  zero  in  the 
formation  of  his  equation  (15),  and  the  usefulness  of  (B.5)  is 
questionable . 

If  the  first  term  on  the  right  hand  side  of  (B.4)  is  replaced 
by  1  and  the  second  and  third  terms  are  divided  by  2,  then  (B.4)  is 
the  appropriate  equation  for  forward  time  differences  instead  of 
centered  time  differences.  If  D  -  0,  this  finite  difference  scheme 
becomes  the  Euler  scheme,  and  it  can  be  shown  that  the  Euler  scheme 
is  unstable  for  almost  all  wavelengths  (Haltiner  and  Williams,  1980, 
p.130).  The  addition  of  diffusion  to  this  method  serves  to 
stabilize  an  otherwise  unstable  numerical  scheme. 

This  stabilizing  effect  is  illustrated  in  Fig.  B.l  where  the 
growth  rate,  G,  is  plotted  against  the  diffusion  parameter,  D.  If  G 
is  greater  than  1  then  numerical  instability  exists.  As  D  increases 
from  0,  G  decreases  until  it  is  less  than  1  and  the  scheme  is 
stable.  As  D  gets  very  large  G  once  again  increases  and  the  method 
becomes  unstable.  So  while  it  is  possible  to  make  the  Euler  scheme 
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stable  by  adding  diffusion,  the  physical  meaning  of  the  diffusion  is 
obscure;  because  the  diffusion  coefficient  must  not  only  represent 
physical  processes,  but  must  also  counteract  the  Inherent  numerical 
instability. 

While  it  is  possible  to  produce  a  finite  difference  formulation 
that  appears  to  be  stable  using  Clancy's  criterion  3.5,  this  is  done 
at  great  computational  expense.  For  Clancy's  ocean  case  (the  case 
presented  in  Fig.  2.1)  Rc  *=  60  and  T  -  .022.  With  F  of  this 
magnitude  it  takes  almost  50  time  steps  for  a  parcel  to  propagate 
through  one  grid  increment.  Using  centered  spatial  differences  the 
time  step  can  be  increased  by  more  than  an  order  of  magnitude  and 
better  results  are  obtained  than  using  the  Clancy  algorithm.  Using 
fourth  order  spatial  differences,  with  a  10  times  larger  time  step, 
the  numerical  results  and  the  analytic  results  are  very  nearly  equal 
(Fig.  B.2).  Any  additional  computations  required  by  the  fourth 
order  scheme  are  more  than  compensated  for  by  the  increased  accuracy 
and  the  larger  time  step. 
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Figure  B.2.  Integration  of  Clancy's  (1981)  oceanic  case  with 

At  *  10  hours.  Same  as  Figure  2.1c  but  with  the 
time  step  increased  an  order  of  magnitude. 

The  numerical  solution  (+)  is  virtually 
indistinguishable  from  the  analytic  solution 
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